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ABSTRACT 


The subject of this research is the buckling behavior of a simply supported 
rectangular plate, with a bulb-flat stiffener attached to one side of the plate. The plate 
structure is subjected to axial compression that increases to the buckling load. The 
stiffener cross-section has a thin web and a bulb-flat flange that extends to one side of the 
web. Results of the investigation include planar property formulas for the asymmetric 
flange geometry, an analytic expression for the Saint Venant torsional constant of the 
flange cross-section, and an analytic expression for the buckling load corresponding to a 
tripping mode of the structure. The torsional constant for the bulb-flat stiffener is 15% - 
23% higher than understood previously. The analytic expression for the buckling load of 
the bulb-flat stiffened plates considered in this investigation yields values that are 2% - 
6% higher than finite element results. It is also shown that the buckling load of a plate 
with a bulb-flat stiffener is 3% - 4% less than that of a plate with a T-fiange stiffener with 
the same cross-sectional area. At the onset of stiffener tripping, the torsionally superior 
bulb-flat tends to bend laterally, while the flexurally superior T-flange tends to twist. 
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EXECUTIVE SUMMARY 


Understanding the elastic stability of stiffened plate structures is of great interest 
in the design of bridges, buildings, automobiles, aircraft, naval vessels, and offshore 
platforms. This work contributes formulas intended for the mechanical engineer, analyst, 
and educator involved in the analysis or design of the bulb-flat cross-section and simply 
supported stiffened rectangular plate. 

Chapter I comments on the intent and synopsis of this investigation. Brief 
comments address the background for the study and literature search for published 
material related to the investigation. A list of notation and conventions is included. 

Chapter II presents the application of multivariable calculus to derive planar 
property formulas for the bulb-flat cross-section. Appendices A and B contain additional 
derivations and expressions related to the formulas. The torsional constant is one of the 
key property values involved in the analysis of stiffened plate structures. No exact 
expression exists in published literature for the torsional constant of the bulb-flat cross- 
section. Saint Venant’s work on an approximate expression for the torsional constant of 
a solid cross-section provides a method to develop an accurate one-term approximate 
formula applicable to bulb-flat cross-sections produced by Corus Group, one of the 
world’s leading producers of bulb-flat profiled metals. Presented is a discussion on the 
idealization method to approximate the torsional constant. 

Chapter III presents the application of elasticity theory and energy methods to 
derive a general expression for the buckling load due to the stiffener tripping of a simply 
supported rectangular stiffened plate subjected to axial compression. The general 
expression for the buckling load involves a constant called . The buckling load and 
values are determined graphically. The value of indicates the stiffener deflection 
behavior at the onset of stiffener tripping. Three special expressions provide initial 
approximations for the upper bound of the critical buckling load value. Concepts applied 
in the derivation include extremum principles in mechanics, the energy criterion, calculus 


XV 



of variations, the Rayleigh-Ritz method, and other topics in applied mathematics and 
physics. The chapter also presents finite element analyses of several stiffened plate 
models. The results are compared to predictions made using derived formulas. MSC 
Nastran 2001/Patran 2001 r3 finite element software and the Maple 8 computer 
environment are the primary analysis and computational tools used in this investigation. 
Appendix C contains data tables used to represent the finite element models. 

Chapter IV summarizes the major findings and conclusions of this work. This 
investigation; 

• Determines boundary equations for the bulb-fiat flange cross-section. 

• Derives planar property value expressions. 

• Determines an approximate torsional constant expression that is more accurate 
than idealizing. 

• Demonstrates that the torsional property of the bulb-flat stiffener is better than 
previously understood. The torque-carrying capacity of a bulb-flat stiffener 
(possessing no structural flaws) is greater than that of an area-equivalent angle 
stiffener. 

• Derives a general expression to predict the buckling load due to stiffener 
tripping of a simply supported rectangular stiffened plate subjected to axial 
compression. The predicted value is less than 6% higher than the finite 
element result. 

• Demonstrates that the buckling behavior of the bulb-fiat stiffened plate is 
unlike that of the T-flange configuration. The bulb-fiat tends to buckle 
laterally and have a lower buckling load value than an area-equivalent T- 
fiange stiffened plate. 

Chapter IV lists future research directions that include: 

• Investigating methods to determine the [X value by other than graphical 
means. 

• Investigating the use of conformal mapping to determine the exact expression 
for the bulb-fiat torsional constant. 

• Investigating solutions to the torsion problem for asymmetric cross-sections. 


XVI 



Conducting investigations of other flange eross-seetions, multiple stiffener 
eonfigurations, and grillages. 

Condueting further investigations of the T flange by treating the flange as a 
thin web plate strip instead of a beam. This treatment may aehieve more 
aeeurate predietions for the T-flange buekling loads. 

Developing algorithms that inerease effieieney in stiffened plate analysis and 
design. 


xvii 
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I. INTRODUCTION 


A. INTENT AND SYNOPSIS 

This work presents an investigation into the struetural behavior and eollapse of a 
stiffened plate panel where the stiffener is thin-webbed and the flange of the stiffener has 
a bulb-flat eross-seetion. Chapter I ineludes brief eomments on the baekground and 
motivation for the study. Chapter II presents analytie expressions for the eross-seetional 
properties of the bulb-flat flange. For brevity, Appendiees A and B eontain key 
derivations related to the property expressions. Cross-seetional properties inelude area, 
eentroid, moments of inertia, and torsional eonstant. The property values for the flange 
eross-seetion are needed in the analytie formulas to determine the buekling load of the 
stiffened plate. Chapter II presents a eomparison of the bulb-flat stiffener eonfiguration 
to the angle flange eonfiguration to show that the torsional property of the bulb-fiat 
stiffener is better than previously understood. Chapter III presents a eomparison of the 
bulb-flat stiffened plate eonfiguration to the T-fiange eonfiguration to show that their 
buekling behaviors are different given the same eross-seetional area. These eomparisons 
provide better diseemment of the buekling eharaeteristies of stiffened plates. Chapter III 
eontains the development of a buekling load formula for a stiffened plate resulting from 
stiffener tripping. The onset of stiffener tripping negates the stiffener’s support to the 
plate panel and leads to eventual eollapse of the strueture. The aim is to develop an 
analytie expression for the buekling load of a reetangular stiffened plate, with one or 
more parallel bulb-flat stiffeners attaehed to one side, where the plate strueture is 
subjeeted to axial eompression. The foeus of the analysis is on a reetangular plate 
stiffened with one bulb-flat stiffener. The results ean be extended and applied to plate 
struetures with more than one stiffener. Energy methods are used, as an alternative to 
veetor methods, to analyze the displaeement and defleetion behavior of the stiffened plate 
model. The ehapter eontains finite element analyses of several stiffened plate models 
using MSC Nastran 2001/Patran 2001-r3 software. Finite element analysis results are 
eompared to formula results to assess formula aeeuraey and the validation of results 
presented earlier. 
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B, BACKGROUND 

Stiffened plates are basie structural components of many items. Metal stiffened 
plates can be found on bridges, buildings, automobiles, aircraft, naval vessels, and 
offshore platforms. In ship design, a stiffener is a metal structure composed of a web and 
flange. Stiffeners placed on one or both sides of a plate add strength and hinder overall 
collapse of the plate panel. Because of the structural complexity of stiffened plates, 
understanding their elastic instability is of great importance to ship designers and requires 
careful study. 

Stiffened plate investigations occur in three general forms: theoretical, using 
classical and emerging theory; numerical, using finite element methods and computer- 
aided simulation; and experimental, using actual grillages. It is enlightening to compare 
the results of all three general forms of investigation. With regard to bulb-fiat stiffened 
plates, it appears that most of what is known of their behavior comes from numerical and 
experimental investigations. This investigation attempts to contribute new theoretical 
insights into the behavior of bulb-flat stiffened plates. 

The use of bulb-flat plate stiffeners in ship design is said to reduce building time 
and maintenance cost. Companies that produce bulb-flat plate stiffeners manufacture 
each stiffener as a single unit, which reportedly results in less production cost compared 
to welded or fabricated stiffeners. Additionally, the curved surface of the bulb-fiat traps 
less moisture resulting in less corrosion. The shape of the bulb-fiat stiffener is much 
easier to inspect, weld, and paint. These benefits save significant repair and maintenance 
costs over the lifetime of the vessel. 

Advocates of bulb-flat structures applaud the advantages gained by using bulb- 
fiats. The Jiangyin Bridge in China and the Oresund Bridge that links Denmark to 
Sweden are major constructions incorporating bulb-fiat geometry. 

C. LITERATURE SEARCH 

In his bibliography, Langhaar includes a list of books and articles that well serve 
the reader interested in general developments and special topics in mechanics. Notable is 
the history and theory of plate stability discussed thoroughly in Bleich [1], who 
acknowledges G.H. Bryan as the originator of the study of plate stability under edge 
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compression. Bleich’s discussion on rectangular plates under uniform axial compression 
and stiffened plates under axial compression are applicable to this investigation. His 
demonstration on the theory of buckling of centrally loaded column s by torsion and 
flexure serves as a fundamental basis for the buckling analysis in Chapter III. 

A series of papers by Chou and Chapman [3] are noteworthy because they address 
the buckling behavior of bulb-flat stiffened plates. Their study presents an improved 
design method, iterative in nature, for determining the buckling load of the structure 
against the torsional buckling of the stiffener. They validate the study using FINASIC, a 
finite element program, and compare the results to actual tests on cruciform struts and 
box columns containing bulb-flat stiffeners. 

An interesting aspect of their theoretical analysis is the idealization of the bulb- 
fiat flange as an equivalent angle flange. That is, they treat the bulb-flat cross-section 
like a rectangular cross-section in regards to the bulb-flat torsional and warping 
properties. Chapter II addresses this treatment and shows that idealization imputes error 
in the calculation of the bulb-flat cross-section’s torsional rigidity. Developing and using 
expressions that maintain the bulb-flat geometry yields more accurate cross-sectional 
property data resulting in a more accurate analysis. 

Other work relevant to this investigation includes research and a series of papers 
by Danielson et al. [5, 6, 7, 8, and 9] analyzing the tripping of a beam attached to a plate 
under lateral pressure loading. His investigation includes the assumption that the 
stiffener behaves like a thin-walled open section beam. Applying nonlinear beam theory, 
he obtains an analytic expression for the buckling of stiffened plate structures under 
longitudinal compression. His work includes investigations of the buckling load of ship 
grillages under axial compression with and without lateral pressure [14 and 15]. He uses 
finite element based eigenvalue analysis to gain insight into the ways the buckling loads 
and modes vary given the grillage dimensions. In new unpublished work, he extends his 
previous work by re-deriving buckling load formulas that incorporate different 
assumptions about buckling behavior. 
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D, 


NOTATION AND CONVENTIONS 


4 

4 

4. 

D 

E 

G 

^ XX ’ ^ yy ’ 4 

4 


4c ’ ^yc 


R 

hf 

K 

K 

m 

r,r^ 

hf 


Area of the flange 
Af+ A , Area of the stiffener 
Area of the web 
Eh^ 


12(1-4) 


, Bending stiffness coeffieient 


Modulus of elastieity of the material (foree/lengm ) 
E 


2(l + v) 


Shear modulus, modulus of rigidity 


Moments of inertia about the axes and eentroid in (x,y) 
coordinate system (length) 

Product of inertia about the origin in {x,y) coordinate system 
(length) 

Moments parallel to the axes through the centroid in {x,y) 
coordinate system (length) 

Product of inertia about the centroid in (x,y) coordinate system 
(length) 

Bending moments (force/length) or first moments (length ) 

Region defined by the cross-section of the flange (length ) 

Height of flange (length) 

Height of the web (length) 

4 + hf, Height of the stiffener (length) 

Slope of a line or line segment 
Radius of curvature (length) 

Thickness of the flange bulb (length) 

4 + 4/, Maximum thickness of flange (length) 
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Thickness of plate (length) 

Thiekness of the web (length) 

Bulb-flat flange eoordinate system 

Bulb-flat flange eoordinate system where the centroid eoineides 
with the origin of the coordinate system 
Plate structure coordinate system 
Angle (radians) 

Poisson’s ratio 

Normal stress (foree/length ) 
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II. BULB-FLAT FLANGE CROSS-SECTION PROPERTIES 


A. INTRODUCTION 

A wealth of knowledge exists on the behavior of stiffeners with common 
geometries such as I, T, Z, angle, channel, and flat-bar cross-sections. These common 
geometries possess certain cross-sectional properties that are relatively easy to derive. 
Studying uncommon cross-sections, such as the bulb-flat cross-section, is often difficult 
because obtaining the section properties of such cross-sections requires considerably 
more effort. 

This chapter presents expressions to determine the planar properties of the bulb- 
flat cross-section. Subsequent analysis of the stiffened plate structure requires knowledge 
of the planar property values for a specified cross-section. 

B, CROSS-SECTIONAL BOUNDARY AND PLANAR PROPERTIES 
Consider the cross-section of a bulb-flat that extends to one side of the web as 

shown in Figures 1-3. Figure 1 shows the entire stiffener cross-section composed of the 
bulb-flat flange and thin web. Figure 2 and Figure 3 show only the flange horizontally 
oriented. Five independent variables {t^, 4^^, r, q, and a ) uniquely determine the 

flange boundary. Using both Figure 2 and Figure 3, the following labeling is established. 
The variable denotes the web thickness defined by the vertical distance between point 
4 and point 5. The variable denotes the bulb thickness defined by the vertical distance 

between point 9 and a point horizontal to point 5. The total flange thickness (at its 
maximum value) is denoted by defined by the vertical distance between point 

9 and the line formed by points 3 and 4. The comers of the cross-section have radii of 
curvature r and q . Points 5-6-7 and 8-9-10 define the corners with radius of curvature 

r. Points 1-2-3 define the corner with radius q. The cross-section has a fiat portion 
defined by the points 7 and 8, that slopes a degrees from the y -axis. 


7 









Figure 2 Bulb-flat flange geometry 



Figure 3 MSC Patran picture of bulb-flat flange boundary 
points for cross-section properties 
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Let denote the flange height. The boundary of the flange ineludes two vertieal sides 


x = 0 and x = hf, where 

/Zy = tan«' +r (l-2tan«' +2see«') (1) 

and Q<a<7r (see Appendix A). The upper and lower boundaries of the eross- 
seetion are defined by the following funetions. The upper boundary of the flange eross- 
seetion is given by: 

fi (^) " (C + hf -r) + ^Jr^ -(x-rf, 0<x<r + rcosa 

/(x) = < /j (x) =-xcotcr + t^+ +r(eot«' + esect-l), r + rcoso!<x<h^ -rcoso! (2) 

/j (x) = + r)-(x-/!y) , hf-rcosa < x< 

The lower boundary of the flange eross-seetion is given by: 


[g2(^) = 0, 


0 < X < Tj 
r^<x<h^ 


( 3 ) 


Given the upper, lower, and vertieal boundaries, the eross-seetion of the flange is 
the bounded region R defined by 0 < x < hj , g(x) <y< /(x), with / (x) and g(x) on 

xG[^0,/2yJ. The flange eross-seetional planar properties ean be ealeulated from the 

following integral expression where m, n are nonnegative integers (see Appendix B for 
details of the integral derivation): 

jj x"/ dA=^ ^ x>" dy dx = —j | x" |[/ (^)]"^ - [g {x)Y } dx (4) 

R 0 g(x) ^ 0 

Substitution of expressions (2) and (3) into expression (4) yields 


^ As CC —> 0 , the flange geometry resembles an angle flange with a rounded end of radius r whose 
eross-seetion properties eould be approximated with an equivalent reetangle. 
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Jjx"= —!— J [r + rsmOy {t^ + t^j-r + rcos,9^ rcosO 

D ^ 1 TT 


de 


hf-rcosa 


where 


If / \«+i 

+- x'” -xcota + r + 4 . + rcota + rcsca-r dx 

n + l J ^ ^ ’ 

r+rcosctr 

1 ^ 

+-j- j (Ay+rsin^) {t^ + r-r cos6)"^\cos6 do 

a-— 

2 

-j (1 + sin (1 - cos cos 9 dO 


hf = hf tan a + r{\-2 tan a+ 2 sec ct) 

0<a<-, --<9<-, 0<x<A 

2 2 2 


/ 


(5) 


Throughout this investigation, the planar property values for the bulb-flat flange are 
caleulated using integral expression (5). Appendix B contains a summary of the planar 
property formulas. For brevity, the expansion of each general planar property formula is 
omitted. 


Integral expression (5) ean be specialized for a elass of bulb-flat flange cross- 
sections by assigning a value to one or more of the variables. For example, by setting 
a = 7t! 6, integral expression (5) defines planar property formulas for a class of bulb-flat 
cross-sections with a 30° slope. Additionally, if the radius of curvature of the eross- 
section eomer is assumed to be one-tenth of the web thiekness (tj =t^/10), the 

integral defines planar property formulas for a specifie type of the 30° slope eross- 
sections. In this case, terms involving beeome negligible for this particular class of 

flange eross-sections. Cross-sectional planar property values for this special class can be 
calculated from the following integral expression (6) where m,n are nonnegative 

integers, a = 7r!6, and the term involving is negleeted: 
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jj x’^y^dA = J [r + rsmOy [t^ + t^^j -r + rcosO'^ rcosOdO 


V3 r, V3 

t^f - \-r 1 h - 

3 6 


+ - 


1 ^ / \/7“l“l 

— I x'" y-xyf3+t^ + ti^j-+ r{\ + yf3^j dx 


n + \ 


1+^ 

2 


1 ° 

+-j- ^ [hf +rsmO^ {t^ + r-r cos0)"^' r cos 6 dO 


(6) 


where 


, V3 

n, =t,f -hr 

j bf 3 


1 + 


2 Vs 


7t 7t 

--<9<-, 0<x</2, 


The following formulas are algebraic expansions of expression (6) for the speeialized 
elass of bulb-flat flanges. The flange area and eentroid formulas are 


4=- 


V " 4y 


2 , 
r + 


1 + 


2 V 3 


L.+ 






Vs 2 Vs 


-bf 


w^'bf 


X = 


n 7t 


Vs 2 


v4 9 9, 


+ 


z 

6 ^ 


r t,., + 


+ 


2 + V 3 


v“ ^ y 
f. rx\ 


^ ^ Vs yrVs^ 

V 


r^t, 


bf 


rtJbf + 


T = 


V y 

^ llVs 7;r 5^ 
^ 18 ~n ~6j 


+ 


1+Vs 

V 6 y 


2 1 2 1 S 

"I- ^w^hf "I- ^hf 

bf ^ w bj 


1 

z 


+ 


V4 “ y 


+ 


5.7r ^ Vs 

I 2 " 




*/ 


y 



^ 2 , 
+ 


rtJbf + 


U"" 3 J 

Vf 

1 3 J 

W 

[2'" 6 J 




Vs 2 Vs 2 Vs 3 

+— t t^,+ — t 4 /+— tJ 

^ w bf ^ W bf 


1 

z 
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Formulas for the moments of the specialized class of flanges are 



I1V3 2;r>/3 2 ^ 

27 9 9 144 


26 Vs 5^ 
27 ^3 


7V3 4 


H— t + 

w of 


9 3 


Vs 3 Vs 4 

+ 27V*/ +jQgV 


i+Vs 


. 2 2 , 

- hf f + 

1 o 


2 V 3 1 


9 3 


43V3 4 2 + V3 

54 3 9 


Vs 1 


Tt r 


H- \tJuf +-1- \tuf K 


27 9 


L-^xx+^yy~ {^^+y^) 

C. TORSIONAL RIGIDITY 

The torsion of a beam is the application of twisting moments at both ends of the 
beam. In the absence of side surface tractions and if the ends of the bar are free to warp, 
the beam experiences pure torsion. The Saint-Venant constant for uniform torsion or 
torsional constant J is the ratio of an applied beam twisting moment to the product 
of the shear modulus (modulus of rigidity) G and the beam twist per unit length 6. 


J = 


GO 


(V) 


The role of the quantity GJ (torsional rigidity) in the twisting of a beam is similar to the 
role of the quantity El (flexural rigidity) in the bending of a beam. 
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1, Exact Value Using the Stress Function Method 

One method of determining the torsional rigidity of a cross-seetion (solving the 
torsion problem) is by finding the stress funetion^. Analysis of torsional behavior 
requires finding the stress funetion O, whieh is an exaet solution to Poisson’s equation: 


^2^ 

v"o = — 


-2G6 in R 


dx^ dy^ 

O = 0 on the boundary of R 


( 8 ) 


Donaldson [10, p. 376] shows the torsional eonstant as 

R 

The following example illustrates using the stress funetion method to determine 
the torsional eonstant. Let the boundary of an elliptie eross-seetion with the origin 
eoineiding with the eentroid be given by 

4+^-i = o (10) 

a b 

whieh has area and inertia properties 

A = 7tab I^=^7lab{a^+b^^ 

Take the stress funetion in the form 


where c is a eonstant. Henee, 


O 


f 2 2 \ 

^ + ^-1 
b^ j 


( 11 ) 


f 


= 2c 


1 1 


\a 


■ 2c 


a^+b^ 

a^b^ 


-2GO 


( 12 ) 


After solving for c, equation (11) satisfies the eonditions of (8) and is zero on the 
boundary of the region when 


2 Introduced by Ludwig Prandtl. 
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Therefore 


c = 


-G9 


+b^ 


J = -- 


V^d) 


jl<l>dA 


2c 


A 

A f 2 2 \ 


dA 


\a 

2 


(1 


2+b2 


' 1 1 

fl 3,^ 

1 1 

fl ,3^ 


—na b 

H—7 

—nab 

-nab 

a^ 

14 ; 


14 7 



2a b 


2u2 


2 , u2 

a +b 


—7tab + —7tab - nab 
4 4 




na b 

2 , r 2 

a +b 


(13) 


When a = b , the results eorrespond with the known formula J = na^ / 2 for the eircular 
eross-seetion. 


The torsion problem has been solved exaetly for other eommon eross-seetional 
forms such as the square, rectangle, and equilateral triangle. The torsional constant can 
be obtained relatively easily for flanges with cross-sections that are symmetric about at 
least one axis. For flanges with arbitrary cross-sections, finding the exact solution is 
quite difficult. The reason for this is, as Donaldson [10, p. 386] notes. 

It is not difficult to write expressions that satisfy the governing differential 
equation, and it is not difficult to write expressions that satisfy the 
boundary conditions. What is difficult is to do both simultaneously. 

The uncommon shape and nonsymmetrical cross-section of the bulb-flat flange make it 

difficult to find a stress function that satisfies simultaneously both conditions of (8). 

Several methods exist that attempt to simplify the search for the stress function of 

“difficult” cross-sections.3 For an arbitrary cross-section, the general solution of the 

torsion problem is found using conformal mapping by mapping the region upon the 

interior and boundary of a circle, then solving the problems of Dirichlet and Neumann for 


3 Methods are generally ineluded in referenees that diseuss solutions to Poisson’s equation or the 
torsion problem. Solving Poisson’s equation using eonformal mapping is diseussed in Henriei [11, pp. 
372 - 377 ]. 
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the circular region. Approximation methods offer alternatives to solving Poisson’s 
equation directly for the bulb-flat flange cross-section. 

2. Approximate Expressions: Lower and Upper Bounds 
Washizu [24] uses variational methods and includes a computational technique 
for obtaining accurate lower and upper bounds for torsional rigidity. His method applies 
the principles of minimum potential and complementary energy. Stakgold [20, pp. 579- 
583] presents a similar approach where he derives 


4 II w dh dy 


IK 


) dx dy 




vj' + ) dx dy 


(14) 


where w vanishes on the boundary of the region and v satisfies 

dx^ dy^ 


(15) 


on the cross-sectional region. The upper bound is smallest if the region’s centroid 
coincides with the origin. Stakgold further shows that the bounds can be improved by 
adding certain harmonic functions to v. It should be noted that the calculations required 
to apply this method to uncommon cross-sections is quite involved. 

The follow examples illustrate establishing bounds using (14) and (15). To 
establish a lower bound for the circular cross-section, let w be the following function 
that vanishes on the circular boundary of radius r whose center coincides with the origin. 

w(x,y) = x^ +y^ (16) 
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A lower bound is 



r 


2 

r 


4 

[[ w dx dy 

4 

II dx dy 


R 

\ ^circle 

) 


V ^circle 

) 


II (w/ + j dx dy II 4x^ + dx dy 




For an upper bound, let 






/ _ 4 



nr 

4 



-Ttr 

) _ 

1 2 

y 


Ttr 


Ttr 


v{x,y) = -Ux^+y^) 


(17) 


that satisfies eondition (15). An upper bound is 

II (v/+v/)jxJj= II (x^ + y^) dx dy = I 

^circle ^circle 

The lower and upper bounds eonfirm the faet that for any eireular eross-seetion, 
I <J < I or J = I that is, the torsional eonstant equals the moment of inertia about the 
eentroid. 


For the square eross-seetion, with sides of length 2a , establish a lower bound by 
ehoosing w to be the following funetion, whieh vanishes on the boundary eentered on 
the origin; 

w{x,y) = {x + a){x-a){y + a){y-a) = [x^-a^) (18) 

A lower bound is 


a a 

4 ^ ^ w dx dy 


Y 

y 




IK 


wj + ) dx dy 


V 9 y 

256a* 

45 


20a' 


:2.2222a' 


(19) 
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An upper bound is established by letting v be the same as equation (17). 

a a a a O 4 

I I [y^ + dx dy = ^ | dx dy = — ^ = 2.6666a‘^ (20) 

-a -a -a -a ^ 

Combining expressions (19) and (20), establishes the following lower and upper bounds 
for the torsional eonstant of a square eross-seetion: 

2.2222a" <J<2.6666a" (21) 

The exaet value of the torsional eonstant for a square eross-seetion is J = 2.2496a'', as 
shown later. 

Applieation of the lower and upper bounds method to the bulb-flat flange eross- 
seetion involves rather arduous ealeulations. Beeause the upper bound is smallest if the 
region’s eentroid eoineides with the origin, the boundary equations (2) and (3) must be 
modified so that the eentroid and origin eoineide. Let the flange eross-seetion dehned by 
equations (2) and (3) be modified for the bounded region (the subseribe c indieating 

the eentroid eoineides with the origin) defined by -x <x^ <hf -x , g{x^) ^ ^ 

with f{x^) and g{x)^ on x^ e [-x,/iy -x] . A funetion w that vanishes on is 

A^c,yc) = [Tc - fl )] [t. - fl (^c )] [Tc - /s (^c )] X 

( 22 ) 

[Tc - (^c )][Tc - ^2 (^c )] (x,+x)(x^-hj-+x) 

Using (22) in (14) results in very eomplieated expressions and provides only modest 
gains to establishing aeeurate bounds for the bulb-flat flange torsional eonstant. Other 
approximation methods may provide simpler means to obtain the torsional eonstant for 
the bulb-flat flange eross-seetion. 
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3, Approximate Expressions: k*A'‘/I 

Saint-Venant^ offers an approximate expression for the torsional constant of any 
solid section (except certain common sections): 




1 

4 ^' 


(23) 


where A is the cross-sectional area and is the moment of inertia about the centroid. 
Let k be defined as 


k = ^ (24) 

A 

for a cross-section where expressions for the area, moment of inertia about the centroid, 
and torsional constant are known. Consider the following table showing the results of 
(24) for common cross-sections. 


Table 1 Saint Venant approximate expression for common cross-sections 


Cross- 

Section 

Parameter 

A 

h 

J 


Circular 

r 

^ 2 

nr 

2 

1 4 

—nr 

2 

^ \ - 0.02533 

4n^r 4n^ 

Elliptic 

a,b 

nab 

^Tlabi^a 

na^b^ 

a^+b^ 

^2 4r4 1 

4 4 4= 2 = 0.02533 

4n a b 4n 

Square 

2a 

4a' 

8 4 

— a 

3 

2.2496a" 

0.02343 


For the given cross-sections, the value of A: is a fixed constant that does not depend on 
the parameters of the cross-section. 

In other cross-sections, the value of k may vary as the parameters that define the 
cross-section vary. For rectangular cross-sections with length 2b and thickness 2t , k 


^ Discussed by Saint-Venant, “Sur une formule donnant approximativement le moment de torsion”, 
Comptes Rendus, vof 88, 1879, pp. 142-154, and in Timoshenko and Goodier [22, pp. 301-302]. 
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varies for different values of the ratio t lb In those cases, the torsional constant can be 

obtained from 



192t 


t ;^tanh 


nTTb^ 



(25) 


developed in Timoshenko and Goodier [22, pp. 309-313]. For a square cross-section, 
(25) yields the J expression shown for the square cross-section in Table 1. For 
rectangular cross-sections with parameters 2t and 2b , the area and inertia properties are 


A = 4bt I^=—bt{b^+t^^ 


The expression for k is thus 

36 




^2 ^ 


1 + - 


V 


V 


1 , riTtb 

. , — tanh- 

^b„=tts,..n 2t 


^ I92t ^ 


Pilkey and Chang [18, p. 103] express the torsional constant formula for rectangular 
cross-sections as 


1 i" 

J--(2h)(2t)'(l-0.63- + 0.052—), where t<b (26) 

which is a simplification of (25). Expression (26) gives adequate accuracy and is used in 
the analysis presented in later sections. For elongated rectangular cross-sections, the 
torsional constant can be obtained from 


J = i(26)(2()= 


where t <Sib 


As the ratio t lb approaches zero, then (26) approaches (27). 


(27) 


^ This effect can be deduced from discussions found in Timoshenko and Goodier [22, pp. 309-313], in 
Donaldson [10, pp.390-394], or in Oden [16, p. 44]. 
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It would be useful to identify elasses or sub-classes of bulb-flat cross-sections 
possessing a fixed constant k . Corns Group^ provides a special profile brochure [4] that 
contains technical data on 59 bulb-flat cross-sections commonly requested. The 
following tables (Table 2 and Table 3) reflect all 59 bulb-flat cross-sections published in 
the brochure. The cross-sections represent a class of bulb-flat flanges where a = 7tl6 
and the values for 4 ^, and r are as given in the tables. The brochure states a 

tolerance level for the value of . For ease of calculations, it is assumed that = C /10, 
which obeys the tolerance level stated in the brochure. 


^ Corus Group was formed in 1999 through the merger of British Steel and Koninklijke Hoogovens. 
The company is a leading international metal company and one of the world’s leading producers of bulb- 
flat profiled metals. 
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Table 2 Finite element results for bulb-flat flange cross-sections 


Variables 

, and r. 

are given, while 

= (i/iox 

and a = n! 6 radians. 

Model 

K 

^bf 

r 

4 

h 

J 

k 


mm 

mm 

mm 

mm^ 

mm^ 

mm" 


1 

6.0 

17.0 

5.0 

335.59 

21903 

13339 

0.023035 

2 

7.0 

17.0 

5.0 

356.15 

24609 

15065 

0.023043 

3 

8.0 

17.0 

5.0 

376.70 

27480 

16915 

0.023084 

4 

6.5 

19.0 

5.5 

410.79 

32907 

19936 

0.023038 

5 

7.0 

19.0 

5.5 

422.19 

34720 

21076 

0.023032 

6 

8.0 

19.0 

5.5 

444.97 

38492 

23475 

0.023049 

7 

10.0 

19.0 

5.5 

490.54 

46667 

28716 

0.023144 

8 

7.0 

22.0 

6.0 

519.50 

53120 

31583 

0.023034 

9 

8.0 

22.0 

6.0 

545.10 

58392 

34823 

0.023031 

10 

9.0 

22.0 

6.0 

570.69 

63911 

38262 

0.023054 

11 

11.5 

22.0 

6.0 

634.66 

78896 

47656 

0.023175 

12 

8.0 

25.0 

7.0 

681.93 

91182 

54641 

0.023039 

13 

9.0 

25.0 

7.0 

711.41 

99088 

59532 

0.023030 

14 

10.0 

25.0 

7.0 

740.89 

107316 

64696 

0.023042 

15 

11.5 

25.0 

7.0 

785.09 

120303 

72912 

0.023088 

16 

8.5 

28.0 

8.0 

849.67 

141249 

85076 

0.023056 

17 

9.0 

28.0 

8.0 

866.36 

146751 

88465 

0.023045 

18 

10.0 

28.0 

8.0 

899.72 

158047 

95505 

0.023035 

19 

11.0 

28.0 

8.0 

933.08 

169751 

102873 

0.023038 

20 

12.0 

28.0 

8.0 

966.43 

181882 

110576 

0.023055 

21 

9.0 

31.0 

9.0 

1035.54 

209455 

126681 

0.023075 

22 

10.0 

31.0 

9.0 

1072.79 

224516 

135976 

0.023049 

23 

11.0 

31.0 

9.0 

1110.03 

240055 

145706 

0.023038 

24 

12.0 

31.0 

9.0 

1147.27 

256095 

155831 

0.023035 

25 

9.5 

34.0 

10.0 

1239.52 

299704 

181896 

0.023095 

26 

10.0 

34.0 

10.0 

1260.02 

309570 

187930 

0.023076 

27 

11.0 

34.0 

10.0 

1301.21 

329720 

200430 

0.023052 

28 

12.0 

34.0 

10.0 

1342.34 

350449 

213450 

0.023039 

29 

10.0 

37.0 

11.0 

1461.61 

416284 

253399 

0.023114 

30 

11.0 

37.0 

11.0 

1506.63 

441894 

269103 

0.023079 
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Table 3 Finite element results for bulb-flat flange cross-sections (cont) 


Variables 

, and r, 

are given, while 

= (1 /10)t^ and a = 7t! 6 radians. 

Model 


V 

r 

4 

h 

J 

k 

31 

12.0 

37.0 

11.0 

1551.64 

468166 

285438 

0.023054 

32 

10.5 

40.0 

12.0 

1701.82 

563855 

344119 

0.023133 

33 

11.0 

40.0 

12.0 

1726.27 

579944 

353915 

0.023113 

34 

12.0 

40.0 

12.0 

1775.17 

612687 

374094 

0.023081 

35 

13.0 

40.0 

12.0 

1824.07 

646209 

395021 

0.023058 

36 

11.0 

43.0 

13.0 

1960.15 

747463 

457230 

0.023151 

37 

12.0 

43.0 

13.0 

2012.93 

787681 

481743 

0.023113 

38 

13.0 

43.0 

13.0 

2065.72 

828771 

507140 

0.023082 

39 

11.5 

46.0 

14.0 

2236.59 

972529 

596102 

0.023167 

40 

12.0 

46.0 

14.0 

2264.93 

997033 

610924 

0.023146 

41 

13.0 

46.0 

14.0 

2321.60 

1046772 

641379 

0.023111 

42 

14.0 

46.0 

14.0 

2378.26 

1097519 

672871 

0.023084 

43 

12.0 

49.0 

15.0 

2531.15 

1244856 

764377 

0.023182 

44 

13.0 

49.0 

15.0 

2591.71 

1304401 

800436 

0.023142 

45 

14.0 

49.0 

15.0 

2652.26 

1365058 

837738 

0.023110 

46 

15.0 

49.0 

15.0 

2712.81 

1426867 

876255 

0.023085 

47 

12.5 

53.5 

16.5 

2990.36 

1736393 

1069001 

0.023213 

48 

13.0 

53.5 

16.5 

3023.55 

1774843 

1092092 

0.023193 

49 

14.0 

53.5 

16.5 

3089.94 

1852703 

1139236 

0.023154 

50 

15.0 

53.5 

16.5 

3156.31 

1931884 

1187887 

0.023123 

51 

16.0 

53.5 

16.5 

3222.69 

2012431 

1237937 

0.023097 

52 

13.0 

58.0 

18.0 

3487.42 

2360220 

1456675 

0.023243 

53 

14.0 

58.0 

18.0 

3559.63 

2458317 

1515332 

0.023202 

54 

15.0 

58.0 

18.0 

3631.84 

2557909 

1575695 

0.023166 

55 

16.0 

58.0 

18.0 

3704.05 

2659049 

1637762 

0.023135 

56 

14.0 

62.5 

19.5 

4061.35 

3199092 

1977224 

0.023249 

57 

15.0 

62.5 

19.5 

4139.39 

3322389 

2050964 

0.023209 

58 

17.0 

62.5 

19.5 

4295.45 

3574247 

2204506 

0.023145 

59 

20.0 

62.5 

19.5 

4529.52 

3966020 

2449749 

0.023082 

] 

Vlean Value 

0.023099 

Standard Deviation 

5.95E-5 

Median Value 

0.023084 
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From Table 2 and Table 3, values for k appear “nearly” fixed for the speeified 
elass of bulb-flat flange eross-seetions. When applied to the bulb-flat flange eross- 
seetion, the Saint Venant approximation (23) with k = \l ~ 0.0253 generates torsion 
values that are 9% to 10% larger than values obtained from finite element methods. 
However, based on the above tables, the following approximate expression yields 
reasonably aeeurate values for the torsional eonstant of flanges within the speeified elass 
when eompared to finite element data. 

A '' 

J,-0.0231^ (28) 

i. 

Henee, expression (28) serves as a working formula that approximates the torsional 
eonstant for the bulb-flange eross-seetions within the speeified elass and is used 
throughout the remainder of this investigation. 

4. Idealization 

Sueeess in determining the exaet expression of Saint Venanf s torsional eonstant 
is limited to simple eross-seetions. It is eommon praetiee to use known exaet expressions 
for simple eross-seetions to approximate more diffieult and multiply eonneeted eross- 
seetions by idealization. Idealizing is attributing a eross-seetion, whose properties may 
be partially known, with the known properties of a simple eross-seetion. Caution is 
required when applying this method to idealize the bulb-fiat flange as an area-equivalent 
angle flange. The following illustrates mathematieally the need for eaution. 

Consider a solid eireular eross-seetion of radius r with area and torsional 
properties, respeetively 


A 

'^circle 


: Ttr J. 


Ttr 


circle 


An area-equivalent square eross-seetion with sides 2/ in length has area and torsional 
properties, respeetively 

4,ua.e=4/^ 4,ua.= 2 . 2496 /^ 
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Since 4.rcle = 4quare ^ then 


/ = 



r 


Hence, 


=2-2496/‘= 2.2496 


'' 

_ 2.2496;r 



8 

V 2 y 


0-8834 (J„.) 


This shows that idealizing a circular cross-section as an area-equivalent square cross- 
section reduces the Saint Venant torsional constant of the cross-section by nearly 12%. 

A similar argument shows that idealizing an elliptical cross-section as an area- 
equivalent rectangular cross-section reduces the Saint Venant torsional constant of the 
cross-section when both cross-sections have high aspect ratios. The aspect ratio for a 
rectangle is the length to thickness ratio Ut. A high aspect ratio means the cross-section 
is very long compared to its thickness. The aspect ratio for an ellipse depends on the 
ratio of the major axis length to the minor axis length. Consider an elliptical cross- 
section centered on the origin with semimajor axis a and semiminor axis b having area 
and torsional properties 


Allipse ‘^ellipse 


Z+F 


An area-equivalent rectangular cross-section of width 21 and thickness 2t , assuming 
t <K /, has area and torsional properties 


A —AlfJ 

'^rectangle rectangle 


\6lt^ 

~Y~ 



Figure 4 Idealizing an ellipse as a rectangle 







As commonly done in this case, assume the two eross-seetions have the same width 
l = a. Sinee , then 

7t 

Alt = Ttab ^ Aat = nab ^t = —b 

4 


Hence, 




^2 y 2 , T 2 ^ 

n a +b 

a" ^ 


na^b^ 

+b^ j 


7t a +b ^ 

< r\ 2 ellipse 

12a j 


The effect on the torsional eonstant due to the idealization depends on the relationship of 
the semimajor and semiminor axis lengths. As the ratio of the semiminor axis to 
semimajor axis approaehes zero, h/a ^ 0 , then 


12 [ a' ^ 


^ — = 0.8225 
12 


Even though methods exist to minimize or eliminate the potential loss of torque eapaeity 
by idealization (detailed discussion omitted), the argument is that application of the 
idealization method requires forethought. 

D, VALIDATION USING FINITE ELEMENT MODELS 

Analysis of several finite element models using MSC Nastran 2001/Patran 2001- 
r3 software serves to validate the property formulas. The property values from formulas 
are compared to values from finite element models. In eaeh finite element bulb-fiat 
flange model, ten distinct points establish the boundary of the bulb-fiat flange eross- 
section. Straight lines and 3-point are curves are defined from the boundary points. A 
single trimmed surfaee is defined from the curves. The arbitrary shape option in the 
software’s beam library is used with the maximum allowable curvature error set at 0.005. 
Maple 8 computer environment and a handheld calculator are used to calculate formula 
results. Table 4 presents planar property values for several bulb-flat flanges using 
expression (28), formula (76) in Appendix A, and formulas (80) - (91) in Appendix B. 


26 



The table presents a comparison of the formula values to data obtained from MSC 
Nastran/Patran finite element models of the cross-sections. 



\ 


\Min Pr 


AREA 

A * 335.5852 
MOMENTS OF INERTIA 

11 =12847.42 

12 = 9055.673 
112 = -3954.662 

TORSIONAL CONST. 
ABOUT CENTROID 
J = 13339.45 
SHEAR STIFF. FACTORS 
K1 = 0.7832503 
K2 = 0.7559787 
WARPING COEFF. 
ABOUT SHEAR CENTER 
CWA = 45319.17 
CENTROID RELATIVE 
TO SHEAR CENTER 
N1A = 0.7564681 
N2A = 0.6324511 
CENTROID RELATIVE 
TO ORIGIN 
HORI. = 8.1627645 
VERT. = 9.3835554 
SHEAR CENTER 
REL.TO ORIGIN 
HORI. = 7.5303135 
VERT. = 8.6270876 
ANGLE FROM AXIS 1 
TO MAX. PRINC. AXIS 
ALPHA = 32.19341 
EXTERNAL PERIMETER 
PERIM = 76.512299 


Figure 5 MSC Patran graph and data of a horizontal hulh-flat flange 
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\Min Pr 


SECTION PROPERTIES 
AREA 

A = 335.586 
MOMENTS OF INERTIA 

11 =15800.79 

12 = 8397.743 

II2 = -3641.825 
TORSIONAL CONST. 
ABOUT CENTROID 
J = 13973.02 
SHEAR STIFF. FACTORS 
K1 =0.6319036 
K2 = 0.1450716 
WARPING COEFF. 
ABOUT SHEAR CENTER 
CWA = 449058.7 
CENTROID RELATIVE 
TO SHEAR CENTER 
N1A = 7.266891 
N2A = 1.493367 
CENTROID RELATIVE 
TO ORIGIN 
HORI. = 7.7303667 
VERT. = 10.266891 
SHEAR CENTER 
REL.TO ORIGIN 
HORI. = 6.2369995 
VERT. = 3. 

ANGLE FROM AXIS 1 
TO MAX. PRINC. AXIS 
ALPHA = 22.26711 
EXTERNAL PERIMETER 
PERIM = 87.175995 

DIMENSIONS 
W= 20.588 
H = 23. 

T1 =6. 

T2 = 12.474 


Figure 6 MSC Patran graph and data of an area-equivalent angle flange 
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Table 4 Comparison of flange values 


Variables , and r, are given, while ry={\l 10)t^ and a = 7t ! 6 radians. 

Percent change = ———xlOO 

* - indicates an absolute percent change of less than 0.01%. 

Source 

K 

^bf 

r 

X 

T 

4 

h 

Percent 

change 

Jf 

Percent 

change 


mm 

mm 

mm 

mm 

mm 

mm- 

mm" 


mm" 


Formula 

6.0 

17.0 

5.0 

8.163 

9.384 

335.59 

21904 

* 

13376 

0.28% 

FE(bulb) 




8.163 

9.384 

335.59 

21903 

* 

13339 

* 

FE(angle) 




7.730 

10.267 

335.59 

24199 

10.48% 

13973 

4.75% 

Formula 

8.0 

17.0 

5.0 

8.397 

10.251 

376.71 

27481 

* 

16928 

0.08% 

FE(bulb) 




8.397 

10.251 

376.70 

27480 

* 

16915 

* 

FE(angle) 




8.010 

11.035 

376.71 

30560 

11.21% 

16026 

-5.26% 

Formula 

IS) 

19.0 

5.5 

9.081 

10.608 

422.19 

34722 

* 

21137 

0.29% 

FE(bulb) 




9.081 

10.607 

422.19 

34720 

* 

21076 

* 

FE(angle) 




8.611 

11.581 

422.19 

38450 

10.74% 

21586 

2.42% 

Formula 

10.0 

19.0 

5.5 

9.408 

11.922 

490.54 

46670 

* 

28661 

-0.19% 

FE(bulb) 




9.408 

11.921 

490.54 

46667 

* 

28716 

* 

FE(angle) 




9.000 

12.754 

490.54 

52004 

11.44% 

26362 

-8.20% 

Formula 

IS 

22.0 

6.0 

10.036 

11.789 

519.51 

53122 

* 

31674 

0.29% 

FE(bulb) 




10.036 

11.789 

519.50 

53120 

* 

31583 

* 

FE(angle) 




9.486 

12.992 

519.51 

58776 

10.65% 

33449 

5.91% 

Formula 

9.0 

22.0 

6.0 

10.287 

12.642 

570.70 

63914 

* 

38339 

0.20% 

FE(bulb) 




10.287 

12.642 

570.69 

63911 

* 

38262 

* 

FE(angle) 




9.784 

13.735 

570.70 

71244 

11.47% 

36965 

-3.39% 

Formula 

8.0 

25.0 

7.0 

11.574 

13.429 

681.94 

91186 

* 

54785 

0.26% 

FE(bulb) 




11.574 

13.428 

681.93 

91182 

* 

54641 

* 

FE(angle) 




10.939 

14.786 

681.94 

100732 

10.47% 

58328 

6.75% 

Formula 

10.0 

25.0 

7.0 

11.829 

14.281 

740.89 

107321 

* 

64856 

0.25% 

FE(bulb) 




11.829 

14.281 

740.89 

107316 

* 

64696 

* 

FE(angle) 




11.243 

15.528 

740.90 

119359 

11.22% 

63541 

-1.79% 


Examination of Table 4 reveals several noteworthy aspeets. The formula results 
agree with the finite element results for the bulb-flat flange. The two independent 
sourees agree almost exaetly on the eentroid loeations (x,y), the flange area values Aj, 

and the values for the polar moments . Though not independent, the formula values for 
the torsional eonstants J^ obtained using the approximate expression (28) differs from 
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the finite element values by less than one half of a percent. Such agreement supports the 
validation of the formula expressions. 

On the other hand, note that the finite element results for the bulb-flat flange 
differ noticeably with the finite element results for an area-equivalent angle flange. The 
centroid coordinates {x,y) differ considerably. The moments differ approximately 

10%. Only the flange area results agree which is by design. This strengthens the case 
that using an angle flange as an area-equivalent to the bulb-flat flange produces error. 

To determine the values for the area and torsional constant for the stiffener that 
includes the web and bulb-fiat flange together, let the stiffener have total height . The 
stiffener area formula is given by 

A,=A,^A, = Aj^(h,-h,)l,. (29) 

The membrane analogy for uniform torsion is a useful tool for visualizing the 
distribution of shearing stress in beam cross-sections and provides justification for 
approximating the torsional constant of the stiffener. If a membrane of constant thickness 
were stretched over a beam cross-section, fixed against the boundary of the cross-section, 
and filled with a gas exerting normal pressure on the membrane surface, the volume 
under the membrane is proportional to the torque-carrying capacity of the cross-section. 
Additionally, the torque-carrying capacity of a cross-section composed of several 
components is greater, but not much greater, than the sum of the components treated 
separately. Hence, the total torsional constant of the entire stiffener is approximated by 
the sum of the component torsional constants (27) and (28) 

A 1 

■f. > Jf +J. = 0.0231-^+-r„= (*. - h,) (30) 

In each finite element stiffener model, eleven distinct points establish the 

boundary of the bulb-flat stiffener cross-section (see Figure 7 and Appendix C). The 

boundary points define straight lines and 3-point arc curves. The curves define a single 

trimmed surface. The beam library provides the cross-sectional properties for the finite 

element models using the arbitrary shape option with the maximum allowable curvature 

error set at 0.005. Table 5 compares the area and torsional constant calculations with 
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finite element model results. Also, Table 5 eompares the area and torsional eonstant 
ealeulations to published data from the Corns teehnieal broehure [4], 




Figure 7 MSC Patran graphs of bulb-flat stiffener boundary points for cross- 
section properties 
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AREA 

A = 932.0546 
MOMENTS OF INERTIA 

11 =1329316. 

12 = 23388.33 
112 = 89132.32 

TORSIONAL CONST. 
ABOUT CENTROID 
J = 20608.63 
SHEAR STIFF. FACTORS 
K1 = 0.6633503 
K2 = 0.1518261 
WARPING COEFF. 
ABOUT SHEAR CENTER 
CWA = 7601410. 
CENTROID RELATIVE 
TO SHEAR CENTER 
N1A = -31.50906 
N2A = 2.364666 
CENTROID RELATIVE 
TO ORIGIN 
HORI. = 5.2984028 
VERT. = 72.076157 
SHEAR CENTER 
REL.TO ORIGIN 
HORI. = 2.9337373 
VERT. = 103.58521 
ANGLE FROM AXIS 1 
TO MAX. PRINC. AXIS 
ALPHA = 176.1135 
EXTERNAL PERIMETER 
PERIM = 275.33551 


Figure 8 MSC Patran graph and data of a bulb-flat stiffener model 


^ax Pr 



SECTION PROPERTIES 
AREA 

A = 932.0578 
MOMENTS OF INERTIA 

11 =1340240. 

12 = 28931.13 
112 = 101280.4 

TORSIONAL CONST. 
ABOUT CENTROID 
J = 18043.08 
SHEAR STIFF. FACTORS 
K1 =0.6610681 
K2 = 0.1392275 
WARPING COEFF. 
ABOUT SHEAR CENTER 
CWA = 8301668. 
CENTROID RELATIVE 
TO SHEAR CENTER 
N1A =-35.13593 
N2A = 2.844582 
CENTROID RELATIVE 
TO ORIGIN 
HORI. = 5.6164331 
VERT. = 72.231934 
SHEAR CENTER 
REL. TO ORIGIN 
HORI. = 2.7718515 
VERT. = 107.36786 
ANGLE FROM AXIS 1 
TO MAX. PRINC. AXIS 
ALPHA = 175.6094 
EXTERNAL PERIMETER 
PERIM = 286. 


Figure 9 MSC Patran graph and data of an area-equivalent angle stiffener model 
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Table 5 Comparison of stiffener values 


Variables h^, t^, 1,^^, and r are given, =(1/10)^^ and a = 7tl6 radians. 

Pereent change = - — xlOO. * - absolute percent change of less than 0.01%. 

Source 

K 

^bf 

r 

K 

4 

"^c-stifffener 

Percent 

change 

*^stiffener 

Percent 

change 


mm 

mm 

mm 

mm 

mm^ 

mm" 


mm" 


Formula 

6.0 

17.0 

5.0 

120.0 

932.06 

1352719 

* 

20534 

-0.37% 

FE(bulb) 





932.05 

1352704 

* 

20609 

* 

FE(angle) 





932.06 

1369171 

1.27% 

18043 

-12.45% 

Coras 





931.00 

1353400 

0.05% 

15950 

-22.61% 

Formula 

8.0 

17.0 

5.0 

120.0 

1172.00 

1676051 

* 

33894 

-0.91% 

FE(bulb) 





1172.00 

1676034 

* 

34204 

* 

FE(angle) 





1172.00 

1694072 

1.08% 

30413 

-11.08% 

Coras 





1170.00 

1681000 

0.30% 

27730 

-18.93% 

Formula 

7.0 

19.0 

5.5 

140.0 

1242.45 

2448920 

* 

34535 

-0.46% 

FE(bulb) 





1242.44 

2448895 

* 

34694 

* 

FE(angle) 





1242.44 

2475783 

1.10% 

30482 

-12.14% 

Coras 





1240.00 

2448000 

-0.04% 

27080 

-21.95% 

Formula 

10.0 

19.0 

5.5 

140.0 

1662.34 

3206519 

* 

67721 

-1.25% 

FE(bulb) 





1662.33 

3206493 

* 

68575 

* 

FE(angle) 





1662.34 

3236434 

0.93% 

61791 

-9.89% 

Coras 





1660.00 

3215600 

0.28% 

57520 

-16.12% 

Formula 

7.0 

22.0 

6.0 

160.0 

1460.10 

3788414 

* 

47037 

-0.26% 

FE(bulb) 





1460.09 

3788380 

* 

47160 

* 

FE(angle) 





1460.10 

3831830 

1.15% 

41325 

-12.37% 

Coras 





1460.00 

3788600 

* 

36810 

-21.95% 

Formula 

9.0 

22.0 

6.0 

160.0 

1780.03 

4558106 

* 

70991 

-0.65% 

FE(bulb) 





1780.02 

4558068 

* 

71453 

* 

FE(angle) 





1780.03 

4605079 

1.03% 

63240 

-11.49% 

Coras 





1780.00 

4553200 

-0.11% 

57630 

-19.35% 

Formula 

8.0 

25.0 

7.0 

180.0 

1885.81 

6188305 

* 

70468 

-0.26% 

FE(bulb) 





1885.80 

6188242 

* 

80674 

* 

FE(angle) 





1885.81 

6261495 

1.18% 

70710 

-12.35% 

Coras 





1890.00 

6189000 

0.01% 

63520 

-21.26% 

Formula 

10.0 

25.0 

7.0 

180.0 

2245.73 

7289798 

* 

115017 

-0.57% 

FE(bulb) 





2245.72 

7289731 

* 

115672 

* 

FE(angle) 





2245.73 

7368371 

1.08% 

102123 

-11.71% 

Coras 





2250.00 

7290500 

0.01% 

93280 

-19.36% 
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In Table 5, each stiffener model contains four rows of results. The four rows are 
referred to as a result set for clarity. Each row is data from a different source. The first 
row in each result set is data from the formulas for the bulb-flat cross-section. The 
second row in each set is data from the finite element analysis for the bulb-flat cross- 
section. The third row in each set is data from the finite element analysis for the area- 
equivalent angle cross-section. And the fourth row in each set is data from the Corns 
technical data brochure. 

Comparing the first two rows of each stiffener result set compares the formula 
results to the finite element results. The independent sources agree almost exactly on the 
stiffener area values and the values for the polar moments /^.^tiffener • The formula 
values for the torsional constants ^stiffener obtained using the approximate expression (30) 
differs from the finite element values by less than two percent. 

Comparing the second and third rows of each stiffener result set compares the 
finite element results of the bulb-flat to the finite element results of the area-equivalent 
angle cross-section. The two sources agree on the stiffener area values and the values for 
the polar moments. However there is a marked difference for the torsional constant 
values. Comparing the second and fourth rows of each stiffener result set compares the 
finite element results of the bulb-flat to the Corus technical data. The sources agree on 
the stiffener area values and the values for the polar moments. The torsional constant 

values differ between 16% and 23%. The reason for the difference appears to be the lack 
of accuracy when using the idealization method to determine the torsional constant value. 

In summary, this chapter presents expressions used to determine the planar 
properties of the bulb-flat cross-section. Stiffened plate analysis requires property value 
accuracy. Five independent parameters, web thickness, bulb thickness, two radii of 
curvature, and the slope angle uniquely define the bulb-flat flange geometry. The 
essential planar property formula (5) takes the form of an integral. All of the property 
values involve integral (5). The conventional methods of determining the torsional 
constant are not viable due to the asymmetrical cross-section of the bulb-flat geometry 
and the complicated equations that define the cross-section boundary. Expression (28) is 
developed from Saint Venanf s one-term approximate expression and serves as a working 
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formula that yields very accurate results compared to finite element methods, thus 
supporting its validation. Using the approximate expression provides values that are 
more accurate than those obtained by idealizing the bulb-fiat flange as an equivalent 
angle flange. The torsional constant for the bulb-flat stiffener is 15% - 23% higher than 
understood previously. In previous studies, the torsional nature of the bulb-flat was 
calculated using the torsional properties of an area-equivalent angle flange stiffener. 
From Table 5, the torque-carrying capacity of the bulb-flat stiffener is found to be greater 
than that of an area-equivalent angle stiffener. 
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III. STIFFENED PLATE BUCKLING 


A, INTRODUCTION 

This chapter deals with the buckling behavior of a stiffened plate. One method of 
analyzing the behavior of struetural members is by deriving the governing differential 
equation and determining the exaet solution using an appropriate differential equation 
teehnique. When the solution is diffieult or impossible to determine, as is the ease with 
stiffened plate struetures, approximate methods of analysis are used. The Rayleigh-Ritz 
method is an effeetive alternate method that yields an approximate solution. The ehapter 
presents an analytie expression for the buckling load of a rectangular stiffened plate, with 
one bulb-flat stiffener attaehed to one side, where the plate structure is subjected to axial 
eompression. The results of the single stiffener analysis can be extended and applied to 
plate structures with more than one bulb-flat stiffener. Additionally, the analytie 
expression ean be used to study the behavior of stiffened plates with different flange 
geometries. 

B, ASSUMPTIONS 


The fundamental assumptions are taken from Danielson [9], modified as 
appropriate, and stated as follows: 


(0 


Eaeh plate-stiffener unit of width b undergoes an identical deformation. 


(ii) 


(iii) 

(iv) 


The plate and web obey the nonlinear Von Karman plate equations (see 
Timoshenko and Gere [21]). The flange obeys the nonlinear beam equations 
derived by Bleieh [1]. 

The plate and stiffener material is elastie, homogeneous, and isotropie. 

Every partiele on the bottom surface of a web undergoes the same 
displaeement as the corresponding particle on the top surface of the plate, and 
every line of partieles in a web normal to the plate surfaee remains normal to 
the deformed plate at its surface. In other words, the bases of the stiffeners 
are clamped to the plate. 
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The prebuckling displacements are less than the maximum thickness of the 
structure and independent of the transverse coordinate. 

The incremental buckling extensional strains at the midsurface of the plate 
and web are negligible. 

The incremental buckling displacements may be approximated by the 
fundamental harmonic in their Fourier expansions. 

The plate and web are so thin that their thicknesses are negligible compared to 
their width, length, and wavelength of deformation. 

The stiffener flange has a solid cross-section with a length that is much greater 
than its largest cross-sectional dimension. Hence, the stiffener flange can be 
treated as a beam. 

The flange undergoes only lateral bending and torsion, no vertical bending. 

C. ENERGY PRINCIPLE 

The energy method for the solution of elastic stability problems is based on an 
extremum principle of mechanics that uses an energy criterion to characterize the 
equilibrium condition of the elastic system. A more precise statement of the energy 
criterion is made in Danielson [5] and is recapitulated in the following discussion. 

Let Pj denote the potential energy of a system in equilibrium state / . Let 

denote the potential energy in a neighborhood of equilibrium state / . Let P\u\ 
represent the increment in potential energy of the elastic system in transition by 
displacement field u from state / to a neighboring state. In addition, let P\u\ be 
expandable into the following component terms (functionals) 

T’[m] = [m] + H2 + ^ "I— 

where Py\u\ refers to linear terms (functionals) with respect to u, refers to 

quadratic terms (functionals) with respect to u , and so forth. Hence 

Pii - Pi = P\u\ = /^ [m] + [m] + /^ [m] H— 


(v) 

(vi) 

(vii) 

(viii) 

(ix) 

(x) 
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If P[u\ > 0 for all non-vanishing neighborhoods of displacement fields u , then 

the fundamental equilibrium state / is stable. In a practical sense, stability here means 
that a force, stress, shock, or disturbance of the system does not cause an excessive 
change or dramatic departure of the equilibrium state to a neighboring state or 
configuration. If there exists a displacement field for which P\u\<0 , the fundamental 

equilibrium state / is unstable. Because the elastic system is in equilibrium, the linear 
term of the potential energy increment must vanish. That is, equilibrium requires 
[m] = 0. The condition necessary for stability is P 2 > 0 • The critical case of neutral 

equilibrium occurs when there exists a displacement field Wj such that [^i ] = 0 and 

^2 [w Mj ] > 0 . In this case, the displacement field Uy is the buckling mode, and the 

value of the load that corresponds to this displacement field is called the bifurcation- 
buckling load. Stated mathematically, 

0<P2[w] (31) 

is the criterion used to determine the buckling mode and load given the properly defined 
total potential energy functional for the stiffened plate structure. 

D, SIMPLY SUPPORTED RECTANGULAR STIFFENED PLATE 

For the stiffened plate structure, the Cartesian coordinate system [xy,X 2 ,x-^) is 
adopted from Danielson’s work to aid result comparison and to minimize confusion with 
the flange cross-sectional coordinate systems {x,y) and Consider a simply 

supported rectangular plate of length a , width b , and thickness , with a longitudinal 

stiffener whose ends are also simply supported. The stiffener divides the width of the 
plate in halves and is composed of a thin web (treated as a plate element) of height 

and thickness , and a bulb-fiat flange (treated as a beam-column element). Suppose the 
stiffened plate structure is under axial compression due to a uniform normal stress <7. 
The quadratic terms of the total potential energy functional for this single stiffened plate 
unit at the instance of buckling can be expressed as 
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( 32 ) 


~ Opiate ^web "*■ ^flange Opiate ^web ^flange 

fV = W[xi,X 2 ) and V = F(xj,X 3 ) are the deflection fields (modes) of the plate 
and stiffener, respectively. Subscripts on W and V denote partial differentiation with 

d^W 

respect to the coordinate system, e.g. =-. In expression (32), U. denotes the 

^x^^X2 

internal potential energy as a result of strains for element i , and T. denotes the potential 

energy of the external loads applied to element i. According to Timoshenko and Gere 
[21, pp. 337, 340, and 350] the internal potential energy expressions due to strains for the 
plate and web are taken as 

b 

1 ^ ^ 

-2{l-v)[w„W^-W„^]}dx,dx,_ (33) 

^ 0 
2 

1 K a 
^ 0 0 

The potential energy expression for the external loads applied to the plate and web are 
taken as 


b 

Opiate KdxA2 (35) 

2 

1 a 

^web = 7 J j V'dx,dx2 (36) 

^ 0 0 

Due to assumption (x), the internal potential energy expression for the flange is 



taken from Timoshenko and Gere [21, p. 25, eq. 1-52 and p. 265, eq. 6-30]. The value of 
I is determined from expression (88). 
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Figure 10 Top view of fiber curvature during buckling 

The potential energy of the external load for the flange is the sum of the products 
of the external forces and the displacements of their points of application in the direction 
of the forces. As the flange buckles, the stresses on the end surfaces may change to 
(7 + d(7. Considering small deformations, the assumption is that the end conditions are 
such that the work done hy d(7 may be neglected in comparison with the work done by 
(7 . Hence the change in potential energy for each fiber in the cross-section of the flange 
is 

= (38) 

Each fiber carries a load crdA and displaces longitudinally a relative distance d^. The 

displacement is due to two reasons: the curvature of the fiber and the change in the 
longitudinal stress. Bleich explains and shows that the potential energy due to the change 
in longitudinal stress vanishes. Hence the change in potential energy for each fiber is due 
to fiber curvature only. The potential energy of the external load for the flange is derived 
here and follows the outline in Bleich [1, pp. 126-127] and Chajes [2, pp. 204-207]. 
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Let denote the relative fiber displaeement due to the eurvature of the flange as 
it buckles under axial compression. The distance is equal to the difference between 
the arc length S of the fiber due to bending and the chord length a of the fiber. Hence 

= S -a 

The potential energy of the external load for the flange is obtained by integrating over the 
cross-sectional area of the flange. Thus 

^flange = CJ-j| S^dx^dx, (39) 

A 

To determine S^, consider a cross-section of the flange at distance x^ along the 
length of the flange. The centroid of the cross-section coincides with the coordinate 
origin. 



Figure 11 Flange cross-section with centroidal coordinate system 

Due to buckling, the point with coordinates {x^,y^) in the cross-section of the flange will 
change coordinates to (x^ + Ax^,y^ + Ay^'j, where Ax^ and Ay^ are functions of Xj. Let 
the point of rotation of the cross-section have coordinates {^,C) relative to the centroid 
of the cross-section. Movements of the point of rotation define the displacement of 
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points in the cross-section: w[x^) in the direction, v(xj) in the direction, and 
/?(xj) counterclockwise rotation about the point of rotation. Hence 


Ax^ cos/? -sin/? x^-^+w 
AyJ |_sin/? cos/? -^ + vj 

Since only small deformations are considered, simplification of the right-hand side gives 

+ w _ x^-^ ^ w-{y^-C)P-vP 
VP 1 j_yc-C+v_ _yc-C_ _v+(x,-^)^+w^_ 

Using the Pythagorean theorem, and considering = Ax^ and = Ay^ to be 

dxi dxy 

small, yields 

^c=‘^-« = n i+(^/) +(^k/) (^0 

The potential energy of the external load for the flange is given by 

W=|crjj| (ax/) +(Ay/) dx,dx^dy^ (40) 

Using the geometrical relations for the area, first and second moments of the flange 
cross-section (assuming the centroid and coordinate origin coincide), 

jj dx,dy^ = 4 

R 

jj x^dx^dy^ = Jj y.dx^dy^ = 0 

R R 

jJ (U - + (U - ^ dx^dy^ = Ij. 

R 

the quadratic terms of the potential energy due to external loads are given by 

^flange (v'Vw'^)-2H^^v'/?'+ 2H^^w'/?' +7^/?'^ Jxi (41) 

A 
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Movements in the direetion affected by >v(xj) are assumed negligible, by assumption 
(x), compared to direction movements perpendicular to the surface of the web 
yielding 

= ^<r] + (42) 

^ 0 

The point of rotation of the cross-section is taken to be the point centered on and at the 
top of the web. This selection provides continuity of the motion at the top of the web and 
the flange. The point of rotation is located half way between points 4 and 5 as shown in 
Figure 3 found on page 9. Relative to the orientation of the flange, ^ is the vertical 
distance from the centroid to the point of rotation. For the bulb-flat cross-section 

^ = hj-X = -{x - hj ^ (43) 

and 



where x is determined from (83) andy is determined from (84). A lateral movement of 
the point of rotation v(xi) equates io V[x^,h^). A rotation or twist about the point of 
rotation of the flange (xj) equates to F 3 (xj, ). 
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1, Energy Functional 

The quadratic terms of the total potential energy functional of a single stiffened 
plate unit at the instance of buckling can be expressed as 


P =77 +77 +77 -T -T -T 

^2 ^ plate ^ web ^ flange plate ^ web flange 


b 

2 a 


11^22 ^12 ^^dx^dX2 


^ 0 
2 


1 

—a 
2 


0 0 
b 

2 a 


11 tpW^dx^dx^ + 11 tj^^dxydx^ +1+ lA^- (x - hj ) V^V^2 +1 ) dx^ 

_b0 0 0 0 * 

2 

The boundary conditions of the stiffened plate structure are 


W{<d,X2) = W{a,X2) = 0 

(46) 

lFii(0,X2) = lFii(a,X2) = 0 

(47) 

F(0,X3) = F (a,X3) = 0 

(48) 

Fii(0,X3) = Fii(a,X3) = 0 

(49) 

F(xi,0) = lF(x„0) 

(50) 

F3(x„0) = lfy(x„0) 

(51) 


2, Approximate Solution 

The Rayleigh-Ritz method is a very effective technique that finds an approximate 
solution to problems requiring the first variation of the total potential energy functional to 
vanish. The method is summarized in two steps. First, assume an admissible solution 
containing unknown coefficients that satisfy the boundary conditions of the problem. 
Second, substitute the assumed solution into the functional and determine the value of the 
unknown coefficients that minimizes the functional. From experimental observation and 
by assumption (vii), the deflection of the plate panel is assumed 
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( 52 ) 


X . mnx, . TtXr, 

XpXj) = csin-i-sin—^ 

a b 

where c is an arbitrary eonstant. The value m is the number of half waves in the 
buekling mode along the longitudinal axis and is taken to be the integer that gives the 
lowest value for the buekling load. The defleetion of the stiffener is assumed 

k (x,, X 3 ) = c sin / (^ 3 ) (53) 

where / (X 3 ) is a funetion to be determined. Expressions (52) and (53) satisfy boundary 

eonditions (46) - (49). Boundary eonditions (50) and (51) transform to the following 
eonditions: 



/( 0 ) = 0 

(54) 

/'( 0 )=f 

(55) 


Henee, 


1 2 2^4 /1 

I c m 7t I 


" 4 ab^ 

Term; =D^^\< 

IT ^ 


/l , 

^ mb a ^ 

( -D„b 

- 1 - 

\2 ^ 

ly a mb J 


+ Term; + f {h^f + GJ^ ^ f {h„f - aJ erm2 

Cl ^ 


t.b^ 


fflTt r / \ Cl / \ 

- f{^3) -/ (^ 3 ) 

a mn 


n2 


+ 2(1-V) /(X3)/^(X3) + /'(X3)" 


\dx. 


Term, =-^ 
' iTl^ 


+ GK\f{^3fdx, + A^f{h^f+2A(x-hf)^f{h^)f{h^) + If^f{h^) 

'77 ' ' 77 77 


h 

n 




(56) 


where the omitted limits of integration are from X 3 = 0 to x-^=h^. 
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Applying the stability criteria (31) yields 


< 


\ ^ mb a ^ 
D b — + — 
2 ^ \ a mb j 


+ Term, + £4 GJ, ^ /'(*, f 


t„b 


l7t^ " 7t 

7t 


+K^\f{x3fdx,+A^.^f{h„f+2Aj\x-h^)^f{h„)f{h„) + If 


Ternii 

TT ^ 


fYlTt r ( \ ^ rff ( \ 

-/(^3)-/ (^ 3 ) 

a mn 


~l2 


+ 2(l-v)l/(x,)/"(x,) + /'(x,) 
(57) 


\dx. 


where a proper selection of / (^ 3 ) minimizes the upper bound for the buckling load. 

Notice that in the case of a plate panel with no stiffener, only the first terms of the 
numerator and denominator in (57) remain and the buckling load value is 


< 


-d/ 
2 ” 


mb a 

— + — 

V a mb 


J _ 


iTT^ 


= D^ 


7t 


2 f 




mb a 
— + — 
V a mb 


Y 

y 


(58) 


where m is one of the two integers closest to alb . Expression (58) agrees with the well- 
known solution for the buckling of a plate that has length a , width b , and is simply 
supported at all edges. 

When there is a thin-webbed stiffener but no flange, the web deflection assumes 
the form 


7 (^ 3 ) 


TtX^ 

T" 


(59) 
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Figure 12 Plate and web deflection when there is no flange 

and the buckling load is the value 



3. General Buckling Mode and Load 

Common conjectures for /(X 3 ) include 


y / \ * 7t%-i 

= —Sin—^ 
h h 

w 

(61) 

r j \ ■ 

/2(x3) = sm 

(62) 


(63) 


Expressions (61) - (63) satisfy boundary conditions (54) and (55). 

,/;(o) = /,(o) = /,(o)=o 
/(o)=/;(o)=/;(o)=| 

When (61) - (63) are substituted into (53), they yield the following stiffener deflection 
approximations: 
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V^{x^,x,) 

V^{x„x,) 

V^{xi,x^) 


= csin 


= csin 


= csin 


mnx^ 

a 

mnx^ 

a 

mnxy 

a 


h . Ttx. 

— sin—i 

b h,., 


sin- 


frx. 


TtX. 


y 


Linear combinations of (61) - (63) are often used to achieve better accuracy in 
approximating the stiffener deflection. The deflection mode function can be assumed in 
the form of a series of functions with undetermined coefficients. If the series of functions 
is complete, by increasing the number of terms in the series the approximation for the 
buckling load converges from above to the exact buckling load value. Functions (61) - 
(63) belong to a family of functions whose members satisfy the boundary conditions. 
Functions belonging to such a family are of the form 


Note that 


/(xj) = ^^^sin^^, 


(64) 


lim/(x 3 )= lim 


kl 


juh^ . nx^ 

—^sm- 




TtX. 


‘w ) 


Continuing with the Rayleigh-Ritz method, substitution of the assumed web 
deflection (64) into the general buckling load expression (57) yields 
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( 65 ) 
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where the value of is chosen to minimize the expression. One method of determining 
/I is by setting equal to zero the ordinary derivative of the quotient expression (65) with 
respect to and solving for . The value can also be obtained graphically. Three 
special values ^ = 1 , = , and ju^°° relate to three special deflection modes. 


For jU = \, the deflection mode becomes f{x^) 


K ■ ^-^3 
— sm—^ 

b K, 


and the buckling 


load is 


< 


-Dpb 

2 


( mb 

a ^ 

' 1 n r 

{ mh. 

a ^ 

— + — 

H— Dh 

+ 

a 

mb ) 

2 w w 

a 

f^K) 


+ GJ, 


In In ^ 


3 


( 66 ) 


Figure 13 depicts the deflection of the structure corresponding to mode 1 that is best 
described by a bending of the web in one half-wave and a rotation of the flange about a 
point at the top of the web. There is no flange bending along its length. Such a mode 
could occur when the flange bending stiffness is large compared to flange torsion and 
web bending. An example of this occurrence is a wide T flange combined with a thin 
web possessing relatively low flexural rigidity. 



Figure 13 Mode 1 deflection 
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For ^ = 2, the deflection mode becomes / (^ 3 ) 


\rl 7tX 

—sin —- and the buckling load is 

b 2h„ 
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(67) 


Figure 14 depicts the deflection of the structure corresponding to mode 2. This deflection 
corresponds to a significant bending of the flange and web with no flange torsion. This 
mode could occur when the flange torsional stiffness is large compared to flange and web 
bending. 



Figure 14 Mode 2 deflection 


/ \ 

For ^ ^ 00 ^ the deflection mode becomes f\x^) = and the buckling load is 
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Figure 15 Mode 3 deflection 

Figure 15 depicts the deflection of the structure corresponding to mode 3. This deflection 
corresponds to the flange exhibiting a combination of bending and twisting while the web 
tends to remain straight. Such a mode could occur when the web flexural stiffness is 
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large compared to flange bending and torsion effects. This case is likely to occur when 
the flange offers little or no additional stiffness to the plate structure compared to the 
contribution of the web. 

The value of [i can be interpreted to indicate the deflection behavior based on the 
three special deflection modes. 

Mode 11 Flange bending Flange torsion Web bending 


1 

1 

No 

Yes 

Yes 

2 

2 

Yes 

No 

Yes 

3 

OO 

Yes 

Yes 

No 
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4, Buckling Expression Summary 

The following summarizes the analytie expressions that predict the critical 
buckling load and buckling behavior of a rectangular stiffened plate with a thin webbed 
stiffener and solid flange treated as a beam. The general expression involves a constant 
fi determined graphically. The value of ii indicates the deflection behavior. 


General Buckling Load Expression 
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Mode 3 Expression (// ^ oo) 
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The three special mode expressions provide initial approximations for the upper bound of 

the buckling load value. Once the proper value of ^ is estimated, the general expression 

provides a value that is less than or equal to the lowest initial approximation. 
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5, Example 

As an example, eonsider a stiffened plate structure, simply supported on all sides, 
subjected to an axial compression load along the longitudinal axis, with the following 
parameters. 

Material properties: 

E = 30,000 ksi v = 0.3 

Plate properties: 

a = 12 inches =1829 mm 
h = 20 inches = 508 mm 
tp =0.3125 inches =7.9375 mm 

Web properties: 

=4.1338 inches =105 mm 
= 0.2362 inches = 6 mm 

Bulb-flat flange properties: 

=0.6693 inches =17 mm 
r = 0.1969 inches =5 mm 
Tj = 0.0236 inches = 0.6 mm 
a=7tl6 radians =30° 


A Maple 8 worksheet provides: 


=83840 
=0.8106 
4=0.0309 
E4 =926071 


4=36212 
4 = 0.5202 
4 =0.0321 
GJf =370796 


m = 4 
x = 0.3214 
4 = 0.2100 


4 = 3.3233 
y = 0.3694 
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The following figure shows the graph of the general expression for the buckling 
load of the plate with a bulb-fiat stiffener. 



Figure 16 Maple 8 graph of the buckling load for the hulh-flat 
stiffened plate example 
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Table 6 Results for the bulb-flat stiffened plate example 


Description 

<7 < Buckling Load 

(psi) 

Remarks 

Plate without stiffener 

26773 

(58) 

Plate stiffened with no flange 

27581 

(60) 

3.02% above the no stiffener 
value. 

Mode 1 

44380 

(66) 

Mode 2 

31262 

(67) 

Mode 3 

31801 

(68) 

General Expression 

27802 

Graphically and (65). 

^ = 3.02 

3.84% above the no stiffener 
value. 

0.80% above the no flange value. 


Without stiffening, the plate is predieted to buckle at a load less than or equal to 26773 
psi. When a flangeless web of height =3.3233 in. is added to the plate panel, the 
plate buckles at a load of less than or equal to 27581 psi. Based on the value for fi , the 
buckling behavior is predicted to exhibit flange bending and torsion with relatively little 
web bending. The given flange cross-section adds relatively little additional stiffness 
compared to the web stiffness. 

E. FINITE ELEMENT ANALYSIS 

Numerical and experimental investigations give insight about the behavior of 
stiffened plates and serve as a basis of comparison to theoretical predictions. The results 
of finite element analyses for several stiffened plate models using MSC Nastran 
2001/Patran 2001-r3 software are compared to analytical predictions to establish a degree 
of confirmation. 
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Figure 17 Stiffened plate pre-buckling and buckling mode 

Several simply supported rectangular plate structures serve as finite element 
models. Each plate structure has a plate panel of length a, width 2b , and thickness 

such that tp«.b<a. The stiffener divides the width of the plate in halves. Each plate 
panel half extends b units of length from the web to the simply supported edge. The 
stiffener includes a thin web of length a, height , and thickness such that 
. Eour flange configurations are examined; no flange, bulb-flat flange, 

circular flange, and T-flange. Each stiffened plate structure is subjected to axial 
compression due to a uniform normal stress < 7 , which increases to the buckling load. 
The boundary of the plate panel is simply supported at all edges. The stiffener is simply 
supported at the ends. 

The derived analytical expressions in the previous section were modeled 
assuming a plate panel of width b . This coincides with the analysis of a single plate unit 
containing a single longitudinal stiffener that would be part of a larger plate structure 
where b units of length separate each stiffener. In this discussion, the plate structure is 
modeled with width 2b , as shown in Eigure 17. The boundary conditions are now 

W (xj,-h) = W (xj,h) = 0 
W22{^v-b) = W22{x„b) = 0 
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Consistency with the finite element model analysis of the 2b width model requires 
changing the limits of integration associated with the plate energy terms and Tpj^jp 

from Xj = ±h / 2 to Xj = ±h . The fundamental problem and its mode shapes remain 
unchanged. 

A secondary aim of the finite element analysis is to estimate the effects of several 
variables or factors. The intention is to discern the sensitivity of the analytic expression 
accuracy due to the variability of web height, flange area, and flange configuration. The 
web height variable is examined at two levels, the flange area at four levels, and the 
flange configuration at four levels. The plate panel parameters for all models remain 
fixed at the following values: 

a = 12 inches =1829 mm 
h = 20 inches = 508 mm 
=0.3125 inches =7.9375 mm 

A 105 mm stiffener with a bulb-fiat flange (see model 1 from Table 2) has a web height 
of 84.41 mm. This is the web height used for the first set of plate models. A 120 mm 
stiffener with a bulb-fiat flange (see model 1 from Table 2) has a web height of 99.41 
mm. This is the web height used for the second set of plate models. The web thickness 
is fixed at 6 mm. 

Web height 1: h^= 84.41 mm 
Web height 2: h^= 99.41 mm 
Web thickness: C = 6 mm 

There are four flange areas examined. The first case is that for a stiffener with no flange. 
The remaining flange area values are derived from the following bulb-fiat flange 
parameters (units in mm) where q = (l/10)c and a = 7t ! 6. 


Flange Area Level 

K 

hf 

r 

Area (mm^) 

1 

N/A 

N/A 

N/A 

0 

2 

6 

17 

5 

335.59 

3 

6 

20 

5 

401.69 

4 

6 

25 

5 

523.39 
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The flange parameters used for area levels 3 and 4 are not from the data in Table 2 or 
Table 3. Analyses of these bulb-flat flanges provide an assessment of the predietion 
accuracy for non-standard parameters. The four flange configurations examined are 


Flange Configuration Label Description 

1 N/A No flange 

2 b Bulb-fiat 

3 c Circular 

4 t T-fiange 

For the stiffened plate structure without a flange, the web deflection mode is of the form 
(59) and the predicted buckling load is determined by 
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Otherwise the stiffener has a flange and the web deflection is of the form (64) with the 
predicted buckling load determined by 
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Using the values from the previous two pages and eonverting the units from 
millimeter to inches provides the parameters for the stiffened plate models. The 
following tables present the converted parameter values. The data tables in Appendix C 
for the stiffened plate finite element models are in inch units (See Table 40 - Table 55). 
Converting the results back into millimeter units may yield values that differ slightly 
from those presented in previous sections. 


Table 7 Summary of stiffened plate model 1 parameters 


Configuration 

Flange Parameters and Properties 

Web Parameters 

Model 1 

4 = 0.5202 in" 

=3.3233 in 

C = 0.2362 in 

No flange 

See Table 40 

No flange area 

Same as above 

Bulb-fiat (b) 

See Table 42^ 

C = 0.2362 in 

= 0.6693 in 

r = 0.1969 in 

=(1/10)C = 0.0236 in 

a = 7t ! 6 radians 
hf = 0.8106 in 

=4.1339 in 

Circular (c) 

See Table 48 

radius = 0.4069 in 
hf = 0.8138 in 

=4.1371 in 

T-flange (t) 

See Table 54 

width = 2.2020 in 

thickness = 0.2362 in 

=3.5595 in 


^ The parameters for this model eonfiguration were used in a grillage tested by the Naval Surfaee 
Warfare Center. The flange and web parameters are the same as used in the example on page 52. 
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Table 8 Summary of stiffened plate model 2 parameters 


Configuration 

Flange Parameters 

Web Parameters 

Model 2 

4 = 0.6226 in" 

= 3.3233 in 

C = 0.2362 in 

No fiange 

See Table 40 

No flange area 

Same as above 

Bulb-flat (b) 

See Table 43 

C = 0.2362 in 

= 0.7874 in 

r = 0.1969 in 

=(1/10X = 0.0236 in 

a = 7t! 6 radians 

= 0.8788 in 

/z, =4.2021 in 

Cireular (e) 

See Table 49 

radius = 0.4452 in 

= 0.8904 in 

=4.2137 in 

T-flange (t) 

See Table 54 

width = 2.6357 in 

thickness = 0.2362 in 

=3.5595 in 


Table 9 Summary of stiffened plate model 3 parameters 


Configuration 

Flange Parameters 

Web Parameters 

Model 3 

4 =0.8113 in" 

=3.3233 in 

C = 0.2362 in 

No flange 

See Table 40 

No flange area 

Same as above 

Bulb-flat (b) 

See Table 44 

C = 0.2362 in 

= 0.9843 in 

r = 0.1969 in 

rj =(1/10)C = 0.0236 in 

a = n! 6 radians 

= 0.9924 in 

\ =4.3157 in 

Circular (c) 

See Table 50 

radius = 0.5082 in 

=1.0163 in 

\ = 4.3396 in 

T-flange (t) 

See Table 54 

width = 3.4343 in = 87.23 mm 

thickness = 0.2362 in = 6 mm 

=3.5595 in 
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Table 10 Summary of stiffened plate model 4 parameters 


Configuration 

Flange Parameters 

Web Parameters 

Model 4 

Same as model 1 

/?„= 3.9138 in 

C = 0.2362 in 

No fiange 

See Table 41 

No flange area 

Same as above 

Bulb-flat (b) 

See Table 45 

Same as model lb 

= 4.7244 in 

Cireular (e) 

See Table 51 

Same as model Ic 

\ = 4.7276 in 

T-flange (t) 

See Table 55 

Same as model It 

= 4.1501 in 


Table 11 Summary of stiffened plate model 5 parameters 


Configuration 

Flange Parameters 

Web Parameters 

Model 5 

Same as model 2 

=3.9138 in 

C = 0.2362 in 

No flange 

See Table 41 

No flange area 

Same as above 

Bulb-flat (b) 

See Table 46 

Same as model 2b 

= 4.7926 in 

Cireular (e) 

See Table 52 

Same as model 2e 

=4.8042 in 

T-flange (t) 

See Table 55 

Same as model 2t 

= 4.1501 in 


Table 12 Summary of stiffened plate model 6 parameters 


Configuration 

Flange Parameters 

Web Parameters 

Model 6 

Same as model 3 

=3.9138 in 

C = 0.2362 in 

No flange 

See Table 41 

No flange area 

Same as above 

Bulb-flat (b) 

See Table 47 

Same as model 3b 

= 4.9063 in 

Cireular (c) 

See Table 53 

Same as model 3e 

=4.9301 in 

T-flange (t) 

See Table 55 

Same as model 3t 

= 4.1501 in 
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The following tables present intermediate results for the stiffened plate models. 


Table 13 Model 1 intermediate results 


Result 

Bulb-flat 

Cireular 

T-flange 

hfiya) 

0.8106 

0.8138 

0.2362 

X (in) 

0.3214 

0.4069 

0.1181 

T(in) 

0.3694 

0.4069 

1.1010 


0.0526 

0.0431 

0.2126 

4 (in') 

0.0309 

0.0215 

0.2102 

4 (in') 

0.0321 

0.0431 

0.0090 

//(in') 

0.2100 

0.1292 

0.2199 

£4(lb-in^) 

926071 

645945 

6305666 

GJ/(lb-m^) 

370796 

496881 

104091 

Mode 1 (psi) 

35599 

36107 

34548 

Mode 2 (psi) 

29053 

28905 

31902 

Mode 3 (psi) 

29356 

29455 

35105 


Table 14 Model 2 intermediate results 


Result 

Bulb-flat 

Circular 

T-flange 

4 (in) 

0.8788 

0.8904 

0.2362 

X (in) 

0.3413 

0.4452 

0.1181 

y(in) 

0.4145 

0.4452 

1.3179 

4 (in') 

0.0773 

0.0617 

0.3633 

/..(in') 

0.0473 

0.0308 

0.3604 

4 (in') 

0.0449 

0.0617 

0.0109 

//(in') 

0.3119 

0.1851 

0.3720 


1417646 

925443 

10813397 

GJ 

517932 

711879 

126078 

Mode 1 (psi) 

36164 

36946 

34614 

Mode 2 (psi) 

29261 

29001 

34229 

Mode 3 (psi) 

30474 

30553 

40794 
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Table 15 Model 3 intermediate results 


Result 

Bulb-flat 

Circular 

T-flange 

hfiya) 

0.9924 

1.0163 

0.2362 

X (in) 

0.3753 

0.5082 

0.1181 

y(in) 

0.4886 

0.5082 

1.7172 

iciyA 

0.1361 

0.1047 

0.8011 


0.0872 

0.0524 

0.7974 


0.0735 

0.1047 

0.0144 

If (in^) 

0.5564 

0.3142 

0.8125 

^4 (lb-in") 

2617311 

1571197 

23921250 

GJ/(lb-in") 

848126 

1208613 

166564 

Mode 1 (psi) 

37428 

38882 

34713 

Mode 2 (psi) 

29798 

29247 

41024 

Mode 3 (psi) 

33128 

33117 

57198 


Table 16 Model 4 intermediate results 


Result 

Bulb-flat 

Circular 

T-flange 

hf{m) 

0.8106 

0.8138 

0.2362 

X (in) 

0.3214 

0.4069 

0.1181 

y(in) 

0.3694 

0.4069 

1.1010 

IcdA 

0.0526 

0.0431 

0.2126 

4(in^) 

0.0309 

0.0215 

0.2102 

Jf{m^) 

0.0321 

0.0431 

0.0090 

/y(in') 

0.2100 

0.1292 

0.2199 

£4 (lb-in") 

926071 

645945 

6305666 

GJ/(lb-in") 

370796 

496881 

104091 

Mode 1 (psi) 

34602 

35109 

33551 

Mode 2 (psi) 

28862 

28657 

32783 

Mode 3 (psi) 

29541 

29495 

37813 
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Table 17 Model 5 intermediate results 


Result 

Bulb-flat 

Cireular 

T-flange 

A/(in) 

0.8788 

0.8904 

0.2362 

X (in) 

0.3413 

0.4452 

0.1181 

y(in) 

0.4145 

0.4452 

1.3179 

4 (in') 

0.0773 

0.0617 

0.3633 

4 (in') 

0.0473 

0.0308 

0.3604 

J,(in') 

0.0449 

0.0617 

0.0109 

If (i4) 

0.3119 

0.1851 

0.3720 

£4(lb-in^) 

1417646 

925443 

10813397 

GJ^(Ib-in^) 

517932 

711879 

126078 

Mode 1 (psi) 

35166 

35947 

33618 

Mode 2 (psi) 

29148 

28791 

35980 

Mode 3 (psi) 

30844 

30672 

45530 


Table 18 Model 6 intermediate results 


Result 

Bulb-flat 

Cireular 

T-flange 

4 (in) 

0.9924 

1.0163 

0.2362 

X (in) 

0.3753 

0.5082 

0.1181 

y(in) 

0.4886 

0.5082 

1.7172 

4 (in') 

0.1361 

0.1047 

0.8011 

4 (in') 

0.0872 

0.0524 

0.7974 

J,(in') 

0.0735 

0.1047 

0.0144 

//(in') 

0.5564 

0.3142 

0.8125 

^/.c(lb-in^) 

2617311 

1571197 

23921250 

GJ^(lb-in^) 

848126 

1208613 

166564 

Mode 1 (psi) 

36430 

37882 

33718 

Mode 2 (psi) 

29886 

29129 

45307 

Mode 3 (psi) 

33964 

33430 

67754 
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The following pictures show the deflection of the plate and web at the onset of 
buckling for Model lb. The plate and web panels buckle in 4 half waves along the 
longitudinal axis. The bulb-flat flange is not visible, but its properties exist in the curve 
that defines the top of the web. 



Figure 18 MSC Nastran picture 1 for stiffened plate Model lb 



Figure 19 MSC Nastran picture 2 for stiffened plate Model lb 
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Circular 

T-flange 


Figure 20 Maple 8 graph 1 of the buckling loads for Model 1 
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Figure 21 Maple 8 graph 2 of the buckling loads for Model 1 
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The pictures below show the difference in the buckling behavior between models 
2b and 3b. In each picture, the top of the web has the beam characteristics of a bulb-flat 
flange. The flange in the left picture has an area property of 0.6226 in . The flange in the 



Figure 23 MSC Nastran pictures showing the buckling behavior of the webs 
with bulb-flat flanges for Models 2b (left) and 3b (right) 

right picture has an area property of 0.8113 in . These properties are not visible in the 
pictures, but exist in the curve that defines the top of each web. The web height in both 
pictures is 3.3233 in. From Table 20, the fi values for the left and right pictures are 2.57 
and 2.25 respectively. This indicates there is more lateral bending of the Model 2b flange 
at the onset of buckling, which agrees with the pictures shown above. 
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The pictures below show the difference in the buckling behavior between models 
4b, 5b, and 6b from left to right respectively. The flanges have area properties of 0.5202 
in , 0.6226 in , and 0.8113 in , respectively. The web height in the pictures is 3.9138 in. 
From Table 20, the // values are 2.68, 2.40, and 2.11, respectively. 



Figure 24 MSC Nastran pictures showing the huckling behavior of 
the wehs with hulh-flat flanges for Models 4h (left), 5h (center), and 6h 
(right) 



Figure 25 MSC Nastran pictures showing 
the huckling behavior of the webs with T-flanges 
for Models 4t (left) and 6t (right) 


The pictures above show the difference in the buckling behavior between models 4t and 
6t. The flanges have area properties of 0.5202 in and 0.8113 in , respectively. The web 
height in the pictures is 3.9138 in. From Table 22, the // values are 1.51 and 1.08, 
respectively. 
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The following tables show the formula predictions and finite element results. 


Table 19 No flange FE results 


Model# 

FE Results (psi) 

Analytic Prediction (psi) 

Percent Change 

1-3 

26928 

27181 

0.94% 

4-6 

26892 

27138 

0.91% 


Table 20 Bulb-flat FE results 


Model# 

FE Results (psi) 

Analytic Prediction (psi) 

A 

Percent Change 

lb 

27672 

28397 

2.86 

2.62% 

2b 

27877 

28824 

2.57 

3.40% 

3b 

28203 

29642 

2.25 

5.10% 

4b 

27663 

28425 

2.68 

2.75% 

5b 

27881 

28908 

2.40 

3.68% 

6b 

28217 

29852 

2.11 

5.79% 


Table 21 Circular FE results 


Model# 

FE Results (psi) 

Analytic Prediction (psi) 

A 

Percent Change 

Ic 

27558 

28246 

2.78 

2.50% 

2c 

27703 

28524 

2.54 

2.96% 

3c 

27936 

29000 

2.29 

3.81% 

4c 

27498 

28178 

2.65 

2.47% 

5c 

27654 

28469 

2.44 

2.95% 

6c 

27907 

28996 

2.21 

3.90% 


Table 22 T-flange FE results 


Model# 

FE Results (psi) 

Analytic Prediction (psi) 

A 

Percent Change 

It 

28767 

31864 

1.85 

10.77% 

2t 

29030 

33200 

1.45 

14.36% 

3t 

29291 

34220 

1.15 

16.83% 

4t 

28739 

32146 

1.51 

11.85% 

5t 

28933 

32945 

1.23 

13.87% 

6t 

29119 

33460 

1.08 

14.91% 
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Table 19 compares the results for the stiffened plates without flanges. The 
formula values nearly agree with the finite element results. Table 20 compares the bulb- 
fiat configurations for each model. The formula values are less than 6% above the finite 
element results. Examination of the fi values indicates that the buckling behavior of the 
stiffened plates tend to have significant bending of the flange and web with a degree of 
flange torsion. This type of buckling behavior can be seen in the Nastran pictures on the 
previous pages. As the flange areas increase between models 1, 2, and 3, and also 
between models 4, 5, and 6, the fj. values tend to approach the value 2, indicating 
perhaps an increasing torsional rigidity effect. The buckling loads tend to increase as the 
flange area increases, as expected. Overall, the results indicate that at the onset of 
buckling, the bulb-fiat stiffened plate structures deflect with a combination of flange 
bending and torsion along with web bending. Table 21 compares the circular flange 
configurations, which appear to buckle at a slightly lower value than the bulb-fiat 
configurations, but behave similar to them at the onset of buckling. Table 22 compares 
the T-flange configurations. From the finite element results, the buckling load of a plate 
with a bulb-fiat stiffener is 3% - 4% less than that of a plate with a T-flange stiffener with 
the same cross-sectional area. Based on the fx values for the T-flange configuration, 
bending of the flange tends to decrease as the cross-sectional area increases. The 
difference in the formula values from the finite element values could be due to 
fundamental assumption (ix). Though the T-flange is solid, treating the flange as a thin 
web plate strip instead of a beam may achieve more accurate predictions. This may 
explain the increase in error for the bulb-fiat and T-flange cases as the area increases due 
to flange widening. The increasing error is less pronounced for the circular flange 
models. 
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IV. FINDINGS AND CONCLUSIONS 


The following summarizes the major findings and eonelusions of this work. Eaeh 
major finding or eonelusion is stated followed by explanatory eomments. 

• Determined eross-seetional boundary equations for the bulb-flat flange eross- 
seetion. 

In order to eonduet an aeeurate analysis of the bulb-flat eross-section, it is 
necessary to determine equations that define or at least approximate the boundary that is a 
closed bounded plane region. The equation that defines the boundary of a simple 
geometric cross-section like a circle or ellipse is well known. The boundary equation for 
uncommon and asymmetric cross-sections is often very difficult to define. As a result of 
this study, functions exist that define the boundary of the bulb-flat flange cross-section. 

• Derived planar property value expressions. 

Determining the critical buckling load of a stiffened plate structure requires 
knowledge of several planar property values of the plate’s cross-section. Based on the 
boundary equations and the application of multivariable calculus, the double integral 
provides expressions for calculating the various planar properties of the bulb-flat flange 
cross-section. Though the integral expressions are complicated, they could be simplified 
by assigning fixed values to certain variables. 

• Determined an approximate torsional constant expression that is more 
accurate than idealizing. 

The Saint-Venant torsional constant is one of the key property values involved in 
the analysis of stiffened plate structures. Due to the uncommon shape and asymmetrical 
property of the bulb-fiat cross-section, determining the exact expression for the Saint- 
Venant torsional constant is difficult and does not exist in published literature. As a 
result of this study, an approximate expression exists for the Saint-Venant torsional 
constant of a specified class of bulb-fiat cross-sections. The approximate expression is a 
one-term function relating the cross-sectional area to the polar moment of inertia. The 
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approximate expression is more accurate than the estimate obtained by idealizing the 
cross-section as an angle flange. 

• Demonstrated that the torsional property of the bulb-flat stiffener is better than 
previously understood. The torque-carrying capacity of a bulb-flat stiffener 
(possessing no structural flaws) is greater than that of an area-equivalent angle 
stiffener. 

In previous investigations the torsional nature of the bulb-flat stiffener was 
idealized as an angle flange stiffener. That is, the bulb-flat flange cross-section was 
treated like a rectangular cross-section in regards to the torsional and warping properties. 
Such treatment imputed error in the calculation of the bulb-flat stiffener’s torsional 
rigidity resulting in conservative estimates. Finite element analysis indicates that the 
torque-carrying capacity of a bulb-flat stiffener (possessing no structural flaws) is greater 
than that of an area-equivalent angle stiffener. 

• Derived a general expression to predict the buckling load due to the stiffener 
tripping of a simply supported rectangular stiffened plate subjected to axial 
compression. In the investigation, the predicted value is less than 6% higher 
than the finite element result. 

Understanding the elastic stability of stiffened plate structures is important to the 
analyst, designer, and educator involved in the analysis and design of structures. Use of 
the energy method provides a technique to derive a general expression for the buckling 
load due to the stiffener tripping of a simply supported rectangular stiffened plate 
subjected to axial compression. As a result of this study, a useful analytic expression 
exists that allows the user to predict the critical buckling load and the buckling behavior 
of a stiffened plate at the onset of stiffener tripping. The onset of stiffener tripping 
negates the stiffener’s support to the plate panel and leads to eventual collapse of the 
structure. The general expression involves a constant called fi that is determined 
graphically. The value of n can be interpreted to indicate the deflection behavior based 
on three special deflection modes. 
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Mode 


Flange bending 

Flange torsion 

Web bending 

1 

1 

No 

Yes 

Yes 

2 

2 

Yes 

No 

Yes 

3 

OO 

Yes 

Yes 

No 


Finite element analyses of several stiffened plate models is presented to validate 
the analytic expressions and assess the sensitivity of the formulas to factor variability. 
Knowledge of MSC Nastran/Patran software and finite element principles along with the 
Maple 8 computer environment was essential for this investigation. 

• Demonstrated that the buckling behavior of the bulb-flat stiffened plate is 
markedly unlike that of the T-fiange configuration. The bulb-fiat tends to 
buckle more laterally and have a lower buckling load value than an area- 
equivalent T-flange stiffened plate. 

The bulb-flat stiffened plate tends to buckle in a Mode 2 fashion, while the T- 
flange stiffened plate tends to buckle in a Mode 1 fashion. Hence, the torsional superior 
bulb-flat tends to bend laterally at the onset of stiffener tripping, while the flexural 
superior T-flange tends to twist axially at the onset of stiffener tripping. 

As a result of this investigation, several essential future research directions exist. 
Future research should include: 

• Investigating methods to determine the [i value by other than graphical 
means. 

• Investigating the use of conformal mapping to determine the exact expression 
for the bulb-flat torsional constant. 

• Investigating solutions to the torsion problem for asymmetric cross-sections. 

• Conducting investigations of other flange cross-sections, multiple stiffener 
configurations, and grillages. 

• Conducting further investigations of the T flange by treating the flange as a 
thin web plate strip instead of a beam. This treatment may achieve more 
accurate predictions for the T-fiange buckling loads. 

• Developing algorithms that increase efficiency in stiffened plate analysis and 
design. 
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APPENDIX A. THE BULB-FLAT FLANGE HEIGHT FORMULA 


This appendix presents the derivation of the formula for the height of the bulb-flat 
flange eross-seetion. Maple 8 eomputer environment software was used to derive some 
of the formula final forms. Five independent variables uniquely determine a speeifie 
bulb-fiat eross-seetion. 

Thiekness of the web (length) 

Thiekness of the flange bulb (length) 

r Radius of eurvature of flange bulb and neek curve (length) 

q Radius of curvature of the cross-section comer (length) 

a Slope angle of the flange neck (radians) 

It will be shown that only three of the independent variables, , r and a , determine the 

height of the flange. Throughout this discussion, it is assumed ^<a<7t H. The 
following figures show a model of the bulb-flat flange oriented horizontally with the bulb 
extended to one side of the web. 
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Figure 26 Bulb-flat flange geometry 



Figure 27 MSC Patran Bulb-flat flange boundary 
points for cross-section properties 
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The three cireles in Figure 26 show the radii of curvature at the corners of the flange. 
The height of the flange, denoted , is the distance from the left boundary x = 0 to the 

right boundary x = h^. Consider the line segment tangent to the two circles of radius r . 

In Figure 27, the line segment connects point 7 ) and point 8 where 


and 


X-, is to be determined 
=C + r-rsina 


Xg = r + r cos a 
y^=K + hf-r + rsma 

The flange height can be expressed simply as 

hj- =x-i+r cos a 

Once Xy is determined, hj- is determined from (71). The slope of line segment is 


(69) 


(70) 


(71) 


, Ay y, - y„ 
slope = — = — -^ = -cot«' 


Ax X 7 - Xg 


•^7 =-(T7 -T8)tana + Xg 


From (69) and (70) 


From (72) and (73) 


and 


y7-y& = {t„ + r-rsma)-[t^ + ti,f-r + rsma) 
= + 2 r(l-sin«') 


X 7 -Xg = tan«'- 2 rtan«' + 2 rsec«'- 2 rcoset 


X 7 = tan«' + r- 2 rtan«' + 2 rsec«'-rcos«' 


Substituting (75) into (71) yields the following formula. 

hy = t^^y tan«' + r (l- 2 tan«' + 2 sec«') 


(72) 

(73) 

(74) 

(75) 

(76) 
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APPENDIX B. THE PLANAR PROPERTY FORMULAS 


The purpose of this appendix is to present the derivation of the planar property 
formulas for the bulb-flat flange. It is assumed ^<a<7t H throughout the discussion. 

Let the closed bounded plane region R be defined by 0 < x < , 

g(x) <y< /(x), with / (x) and g(x) on x G [O, hj ] , such that 


and 


where 


f{^) 


f, (x) = , 0<x<r + rcosa 

f^[x) = mx + b, r + rco^a<x<h^-rcosa 


fi ix)^(t„+r)-^r"-(x-hf )', hj.-r 


cosa<x<hj 


g{x) 


g 2 (x) = 0, Tj < X < hj. 


< r. 


m = -cotcr 

b = t^ + r i^coia + csca - \) 
hf =t^f tan«'-l-r(l-2tan«'-l-2sec«') 


(77) 


(78) 


(79) 


By the theorem of integrability of a piecewise function*, since / and g are piecewise 
defined functions on the same closed interval, then / and g are integrable on the same 
closed interval. The planar property formulas can be determined from the following 
expression where m,n are nonnegative integers. 

jjx”’y'’dA=j j x>VyJx = ^--Jx'”|[/(x)]"^‘-[g(x)]"^'| Jx 

« 0 g(x) ^ 0 ^ 


* Ross, K.A., Elementary Analysis: The Theory of Calculus, Springer, New York, pp. 258-259. 
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—-cir 

jjxy"ciA = ^ j (r + rsinffj'^ (^t^+t^j^-r + rcosffj rcosOdO 

R ^ + 1 

2 

hf-rcosa 

If ( \«+l 

+- x'” I -X cot ct ++ 4 . + r cot ct + r CSC ct - r) dx 

n + 1 ^ ' 

/■+/■ coster 

1 ^ 

+-j- J (/Zy +rsin^) {t^ + r-r cosr cos6 dO 


-'i 

n + \ 


I {\ + smOy {l-cosO)"^'^ cosO do 


where 

hf = hf tan ct + r (l - 2 tan ct + 2 sec a) 

0 <ct<-, --<6<-, 0<x<hi 

2 2 2 ^ 

By definition, the area of the closed bounded plane region R is the value of the integral 

hf f(x) hf 

Area = || dA = ^^dA = ^ J dy dx = ^ f i^x)- gi^x) dx 

R R 0 g{x) 0 

Hence, 

r+rcoso: hf-rcosa hf hf 

II dH = I /j(x)d'x+ I f2{x)dx+ I /j (x)d'x-|gj (x)d'x-I gj 


2 / ~ 

jjdA= I +tj,^-r+ reos6')reos^ | (-xcotcr + t^+ 4^+rcot«' + rcse«'-r)d'x 


+ 1 1... + r-r 


eos6')r eos6' dO-r^ | (l-eos^)cos^ dO 


Integrating and collecting terms yields the area of the bulb-flat flange: 


Af=- 1-^ (1 - 2 tancr + 2sec«') +(1 - tancr + secct) 


1 2 

H— t., tan a + tj,,. tan a 


82 



Other planar property formulas are determined in similar fashion and are 
summarized below (algebraic expansion omitted). 


= jj x^y^dA = ^^ydA 

R R 

n 

—a 

1 ^ 


cos^) rcosOdO 


+— I (-xcotcr + t^+ +rcotcr + rcsccr-r) dx 

^ r+rcosdtr 

1 0 ^0 

+ — j [t^ + r-rcosO)^ rcosO do — | (l-cos^)^ cos6^ 

n ^ n 

a — — 

2 2 

“ II “ II ^ 

R R 

n 

—a 
2 

= j (r+ rsm^)(t^+ -r+ rcos^)rcos6'^6^ 


j xi^-xcoia + t^+tf^j- +rcot«' + rcsc«'-r) dx 


+ rsmO' (+ r-r 


u 

cos r cos 9 d9-r^ ^ (l + sin (l - cos cos 9 d9 


83 



= JJ x^y^dA = Jj dA 

R R 

n 

—a 

1 / >.3 

= - J [t^+tf^f -r + rcosd^ rcos 0 d 0 


+ - j {^-xcota + t^ + ti^f^ +rcota + rcsca-r^ dx 

^ r+rcosctr 

1 0 .0 
+ ^ I -rcos^)^ r cos^ d0 — j (l-cos^)^ cos0 d0 

n n 

a — — 

2 2 

lyy = jJ x^y^dA = Jj dA 

R R 

71 

—a 

2 ^ 

= j (r + rsin^) -r+ rcos0^rcos0 d0 

71 

~i ^ 

hf -rcosa 

+ J x^ {-xcoia + t^+tf^j^+rcoia + rcsca-r^dx 

r+rcosctr 

0 2 ^ 

+ j [hj +rsm0) [t^ + r - r COS0) r cos 0 d0-r^ j (l + sin^)^ (l -cos0)cos0 d0 

71 71 

a — — 

2 2 

" jj = jj dA 

R R 

71 

--or 

= — J {r + rsm0^[t^+tj^j-r + rcos0^ rcos0 d0 

^ h^-rco^a 

+ — J x(-xcot«' + ^^+rcot«' + rcsc«'-r) dx ( 

r-\-rcosa 

1 ^ 2 

+ 9 j (^/ + ''sin0)(^„ + r -rcos0) rcos0d0 

2 JZ 

a — 

2 

1 ^ 

— Tj"* j (l + sin6')(l -0086^)^ COS0 d0 
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( 89 ) 


^xyc=hy-My ( 90 ) 

4=4+^.c (91) 
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APPENDIX C. DATA TABLES 


This appendix presents flange, stiffener, and plate data tables used to model the 
cross-seetions for finite element analysis. When creating a finite element model of a 
structure, a matrix represents the structure. The information in the left portion of each 
table list parameters or references related to the set of points contained in the right 
portion of the table. 

1, DATA FOR BULB-FLAT FLANGE CROSS-SECTION 

See Figure 3 and Figure 5. 


Table 23 Data for Flange Model 1 in Table 2 (units in mm) 


= 6 

[0.0000 0.6000 0.0000][0.1757 0.1757 0.0000][0.6000 0.0000 

=17 

0.0000][20.5885 0.0000 0.0000][20.5885 6.0000 0.0000][18.0885 


6.6699 0.0000][16.2583 8.5000 0.0000][9.3301 20.5000 

r = 5 

0.0000][5.0000 23.0000 0.0000][0.0000 18.0000 0.0000] 

II 

o 

II 

ON 

o 


a = 7r 16 


hf. = 20.59 


4 =335.59 



Table 24 Data for Flange Model 3 in Table 2 (units in mm) 


t =8 

w 

[0.0000 0.8000 0.0000][0.2343 0.2343 0.0000][0.8000 0.0000 

6. =17 

0.0000][20.5885 0.0000 0.0000][20.5885 8.0000 0.0000][18.0885 

bf ^ ' 

8.6699 0.0000][16.2583 10.5000 0.0000][9.3301 22.5000 

r = 5 

0.0000][5.0000 25.0000 0.0000][0.0000 20.0000 0.0000] 

o 

So 

II 

II 


a = 7t/6 


hf = 20.59 


Af =316.11 
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Table 25 Data for Flange Model 5 in Table 2 (units in mm) 



[0.0000 0.7000 0.0000][0.2050 0.2050 0.0000][0.7000 0.0000 

t.r =19 

0.0000][22.8205 0.0000 0.0000][22.8205 7.0000 0.0000][20.0705 


7.7369 0.0000][18.0574 9.7500 0.0000][10.2631 23.2500 

r = 5,5 

0.0000][5.5000 26.0000 0.0000][0.0000 20.5000 0.0000] 

II 

o 

II 

o 


a = 7t/6 


= 22.82 


4 =422.19 



Table 26 Data for Flange Model 7 in Table 2 (units in mm) 


^w=10 

[0.0000 1.0000 0.0000][0.2929 0.2929 0.0000][1.0000 0.0000 

6. =19 

0.0000][22.8205 0.0000 0.0000][22.8205 10.0000 

bf 

0.0000][20.0705 10.7369 0.0000][18.0574 12.7500 

r = 5,5 

0.0000][10.2631 26.2500 0.0000][5.5000 29.0000 

r,=(l/10)^„=10/10=l 

0.0000][0.0000 23.5000 0.0000] 

a = 7t/6 


hf = 22.82 


4 = 490.54 



Table 27 Data for Flange Model 8 in Table 2 (units in mm) 



[0.0000 0.7000 0.0000][0.2050 0.2050 0.0000][0.7000 0.0000 

tr = 22 

0.0000][25.6299 0.0000 0.0000][25.6299 7.0000 


0.0000][22.6299 7.8038 0.0000][20.4338 10.0000 

r = 6 

0.0000][11.1962 26.0000 0.0000][6.0000 29.0000 

II 

o 

II 

o 

0.0000][0.0000 23.0000 0.0000] 

a = 7tl6 


hf = 25.63 


4 =519.51 
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Table 28 Data for Flange Model 10 in Table 2 (units in mm) 


K=9 

[0.0000 0.9000 0.0000][0.2636 0.2636 0.0000][0.9000 0.0000 

= 22 

0.0000][25.6299 0.0000 0.0000][25.6299 9.0000 


0.0000][22.6299 9.8038 0.0000][20.4338 12.0000 

r = 6 

0.0000][11.1962 28.0000 0.0000][6.0000 31.0000 

II 

o 

II 

o 

0.0000][0.0000 25.0000 0.0000] 

a = 7t/6 


hf = 25.63 


Af = 570.70 



Table 29 Data for Flange Model 12 in Table 2 (units in mm) 


t =8 

w 

[0.0000 0.8000 0.0000][0.2343 0.2343 0.0000][0.8000 0.0000 

6. = 25 

0.0000][29.5167 0.0000 0.0000][29.5167 8.0000 

bf 

0.0000][26.0167 8.9378 0.0000][23.4545 11.5000 

r = 1 

0.0000][13.0622 29.5000 0.0000][7.0000 33.0000 

o 

So 

II 

II 

0.0000][0.0000 26.0000 0.0000] 

a = 7tl6 


hf = 29.52 


4 = 681.94 



Table 30 Data for Flange Model 14 in Table 2 (units in mm) 


^w=10 

[0.0000 1.0000 0.0000][0.2929 0.2929 0.0000][1.0000 0.0000 

6. = 25 

0.0000][29.5167 0.0000 0.0000][29.5167 10.0000 


0.0000][26.0167 10.9378 0.0000][23.4545 13.5000 

r = l 

0.0000][13.0622 31.5000 0.0000][7.0000 35.0000 

r,=(l/10)^„=10/10=l 

0.0000][0.0000 28.0000 0.0000] 

a = 7tl6 


hf = 29.52 


4 = 740.89 
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2 . 


DATA FOR THE BULB-FLAT STIFFENER CROSS-SECTION 


See Figure 1, Figure 7, and Figure 8. 


Table 31 Data for first stiffener in Table 5 (units in mm) 


o 

(N 

II 

[0.0000 0.0000 0.0000][6.0000 0.0000 0.0000][6.0000 99.4115 

hr = 20 59 

0.0000][6.6699 101.9115 0.0000][8.5000 103.7417 


0.0000][20.5000 110.6699 0.0000][23.0000 115.0000 

99.41 

0.0000][18.0000 120.0000 0.0000][0.6000 120.0000 

^w = 6 

0.0000][0.1757 119.8243 0.0000][0.0000 119.4000 0.0000] 



r = 5 


II 

o 

II 

ON 

o 


a = 7tl6 



Table 32 Data for second stiffener in Table 5 (units in mm) 


A, =120 

[0.0000 0.0000 0.0000][8.0000 0.0000 0.0000][8.0000 99.4115 

hr =20 59 

0.0000][8.6699 101.9115 0.0000][10.5000 103.7417 


0.0000][22.5000 110.6699 0.0000][25.0000 115.0000 

99.41 

0.0000][20.0000 120.0000 0.0000][0.8000 120.0000 


0.0000][0.2343 119.7657 0.0000][0.0000 119.2000 0.0000] 



r = 5 


o 

So 

II 

II 


a = 7r 16 
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Table 33 Data for third stiffener in Table 5 (units in mm) 


= 140 

[0.0000 0.0000 0.0000][7.0000 0.0000 0.0000][7.0000 117.1795 

hf. = 22.82 

0.0000][7.7369 119.9295 0.0000][9.7500 121.9426 


0.0000][23.2500 129.7369 0.0000][26.0000 134.5000 

=117.18 

0.0000][20.5000 140.0000 0.0000][0.7000 140.0000 


0.0000][0.2050 139.7950 0.0000][0.0000 139.3000 0.0000] 



r = 5.5 


II 

o 

II 

o 


a = 7t/6 



Table 34 Data for fourth stiffener in Table 5 (units in mm) 


= 140 

[0.0000 0.0000 0.0000][10.0000 0.0000 0.0000][10.0000 

h = 22 82 

117.1795 0.0000][10.7369 119.9295 0.0000][12.7500 121.9426 


0.0000][26.2500 129.7369 0.0000][29.0000 134.5000 

=117.18 

0.0000][23.5000 140.0000 0.0000][1.0000 140.0000 

^.=10 

0.0000][0.2929 139.7071 0.0000][0.0000 139.0000 0.0000] 

V=19 


r = 5.5 


r,=(l/10)^„=10/10=l 


a = 7t/6 



Table 35 Data for fifth stiffener iuTable 5 (units in mm) 


\ = 160 

[0.0000 0.0000 0.0000][7.0000 0.0000 0.0000][7.0000 134.3701 

h, = 25 63 

0.0000][7.8038 137.3701 0.0000][10.0000 139.5662 


0.0000][26.0000 148.8038 0.0000][29.0000 154.0000 

A„= 134.37 

0.0000][23.0000 160.0000 0.0000][0.7000 160.0000 


0.0000][0.2050 159.7950 0.0000][0.0000 159.3000 0.0000] 

hf ~ 


r = 6 


II 

o 

II 

o 


a = 7t/6 
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Table 36 Data for sixth stiffener in Table 5 (units in mm) 


= 160 

[0.0000 0.0000 0.0000][9.0000 0.0000 0.0000][9.0000 134.3701 

h. =25 63 

0.0000][9.8038 137.3701 0.0000][12.0000 139.5662 


0.0000][28.0000 148.8038 0.0000][31.0000 154.0000 

134.37 

0.0000][25.0000 160.0000 0.0000][0.9000 160.0000 

^w=9 

0.0000][0.2636 159.7364 0.0000][0.0000 159.1000 0.0000] 

(N 

(N 

II 


r = 6 


II 

o 

II 

o 


a = 7tl6 



Table 37 Data for seventh stiffener in Table 5 (units in mm) 


o 

oo 

II 

[0.0000 0.0000 0.0000][8.0000 0.0000 0.0000][8.0000 150.4833 

hf = 29.52 

0.0000][8.9378 153.9833 0.0000][11.5000 156.5455 


0.0000][29.5000 166.9378 0.0000][33.0000 173.0000 

150.48 

0.0000][26.0000 180.0000 0.0000][0.8000 180.0000 


0.0000][0.2343 179.7657 0.0000][0.0000 179.2000 0.0000] 

(N 

II 


r = l 


o 

oo 

II 

II 


a = 7tl6 



Table 38 Data for eighth stiffener in Table 5 (units in mm) 


o 

oo 

II 

[0.0000 0.0000 0.0000][10.0000 0.0000 0.0000][10.0000 

hf = 29.52 

150.4833 0.0000][10.9378 153.9833 0.0000][13.5000 156.5455 


0.0000][31.5000 166.9378 0.0000][35.0000 173.0000 

150.48 

0.0000][28.0000 180.0000 0.0000][1.0000 180.0000 

^w=10 

0.0000][0.2929 179.7071 0.0000][0.0000 179.0000 0.0000] 

hf ~ 25 


r = l 


r,=(l/10)^„=10/10=l 


a = 7tl6 
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Table 39 Data for angle stiffener in Figure 9 (units in mm) 


A, =120 

[0 0 0][6 0 0][6 107.5259068 0][23 107.5259068 0][23 120 0][0 120 0] 

^..=107.53 


C = 6 




(N 

II 



3. DATA FOR THE STIFFENED PLATE MODELS 


Table 40 Data for stiffened plate models 1-3 with no flange (units in inches) 


See Table 7 for 
parameters. 

[0.0000 0.0000 3.3233][72.0000 0.0000 3.3233][0.0000 0.0000 
0.0000][72.0000 0.0000 0.0000][0.0000 -20.0000 0.0000][72.0000 
-20.0000 0.0000][0.0000 20.0000 0.0000][72.0000 20.0000 0.0000] 

Table 41 Data for stiffened plate models 4-6 with no flange (units in inches) 

See Table 10 for 
parameters. 

[0.0000 0.0000 3.9138][72.0000 0.0000 3.9138][0.0000 0.0000 
0.0000][72.0000 0.0000 0.0000][0.0000 -20.0000 0.0000][72.0000 
-20.0000 0.0000][0.0000 20.0000 0.0000][72.0000 20.0000 0.0000] 


Use the following tables to build the stiffened plate models. In MSC Patran, 
create the geometry from the points in each table by entering the points and then create 
the curves and surfaces as appropriate. Create a mesh (with mesh size 1) for the plate and 
web surfaces. Create a mesh for the curve at the top of the web with topology BAR2. 
The curve will serve as the flange once the appropriate beam section is created. Using 
the properties menu, create and name a ID beam property set. Use the general section 
(CBEAM) option. Use the input properties button and then the create sections beam 
library to create and emplace the appropriate flange onto the web. For the bulb-flat 
flange, use the beam library to create and name a bulb-flat flange by selecting the 
arbitrary shape button and the boundary loops method button. For the input option, set 
the select surface option. In the select surface menu, set the maximum allowable 
curvature error to 0.005. Select the surface and press OK. In the general beam 
(CBFAM) input properties subordinate menu, enter the appropriate material name. 
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section name (just created in the create sections beam library), and bar orientation (start 
by using the vector (0,1,0) and adjust with (l,0,0) or (0,0,l) as needed). 


Table 42 Data for bulb-flat stiffened plate model lb (units in inches) 


See Table 7 model 
lb for parameters. 

[0.0000 -0.1181 3.3233][0.0000 0.1181 3.3233][0.0000 0.1445 
3.4217][0.0000 0.2165 3.4938][0.0000 0.6890 3.7665][0.0000 

0.7874 3.9370][0.0000 0.5906 4.1339][0.0000 -0.0945 
4.1339][0.0000 -0.1112 4.1269][0.0000 -0.1181 4.1102][0.0000 
0.0000 3.3233][72.0000 0.0000 3.3233][0.0000 0.0000 
0.0000][72.0000 0.0000 0.0000][0.0000 -20.0000 0.0000][72.0000 
-20.0000 0.0000][0.0000 20.0000 0.0000][72.0000 20.0000 0.0000] 

Table 43 Data for bulb-flat stiffened plate model 2b (units in inches) 

See Table 8 model 
2b for parameters. 

[0.0000 -0.1181 3.3233][0.0000 0.1181 3.3233][0.0000 0.1445 
3.4217][0.0000 0.2165 3.4938][0.0000 0.8071 3.8347][0.0000 

0.9055 4.0052][0.0000 0.7087 4.2020][0.0000 -0.0945 
4.2020][0.0000 -0.1112 4.1951][0.0000 -0.1181 4.1784][0.0000 
0.0000 3.3233][72.0000 0.0000 3.3233][0.0000 0.0000 
0.0000][72.0000 0.0000 0.0000][0.0000 -20.0000 0.0000][72.0000 
-20.0000 0.00001[0.0000 20.0000 0.00001[72.0000 20.0000 0.0000] 

Table 44 Data for bulb-flat stiffened plate model 3b (units in inches) 

See Table 9 model 
3b for parameters. 

[0.0000 -0.1181 3.3233][0.0000 0.1181 3.3233][0.0000 0.1445 
3.4217][0.0000 0.2165 3.4938][0.0000 1.0039 3.9484][0.0000 

1.1024 4.1189][0.0000 0.9055 4.3157][0.0000 -0.0945 
4.3157][0.0000 -0.1112 4.3088][0.0000 -0.1181 4.2921][0.0000 
0.0000 3.3233][72.0000 0.0000 3.3233][0.0000 0.0000 
0.0000][72.0000 0.0000 0.0000][0.0000 -20.0000 0.0000][72.0000 
-20.0000 0.0000][0.0000 20.0000 0.0000][72.0000 20.0000 0.0000] 

Table 45 Data for bulb-flat stiffened plate model 4b (units in inches) 

See Table 10 
model 4b for 
parameters. 

[0.0000 -0.1181 3.9138][0.0000 0.1181 3.9138][0.0000 0.1445 
4.0123][0.0000 0.2165 4.0843][0.0000 0.6890 4.3571][0.0000 

0.7874 4.5276][0.0000 0.5906 4.7244][0.0000 -0.0945 
4.7244][0.0000 -0.1112 4.7175][0.0000 -0.1181 4.7008][0.0000 
0.0000 3.9138][72.0000 0.0000 3.9138][0.0000 0.0000 
0.0000][72.0000 0.0000 0.0000][0.0000 -20.0000 0.0000][72.0000 
-20.0000 0.0000][0.0000 20.0000 0.0000][72.0000 20.0000 0.0000] 
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Table 46 Data for bulb-flat stiffened plate model 5b (units in inches) 

See Table 11 [0.0000-0.1181 3.9138][0.0000 0.1181 3.9138][0.0000 0.1445 

model 5b for 4.0123][0.0000 0.2165 4.0843][0.0000 0.8071 4.4253][0.0000 

parameters. 0.9055 4.5958][0.0000 0.7087 4.7926][0.0000 -0.0945 

4.7926][0.0000 -0.1112 4.7857][0.0000 -0.1181 4.7690][0.0000 
0.0000 3.9138][72.0000 0.0000 3.9138][0.0000 0.0000 
0.0000][72.0000 0.0000 0.0000][0.0000 -20.0000 0.0000][72.0000 
_ -20.0000 0.00001[0.0000 20.0000 0.00001[72.0000 20.0000 0.0000] 


Table 47 Data for bulb-flat stiffened plate model 6b (units in inches) 

See Table 12 [0.0000-0.1181 3.9138][0.0000 0.1181 3.9138][0.0000 0.1445 

model 6b for 4.0123][0.0000 0.2165 4.0843][0.0000 1.0039 4.5389][0.0000 

parameters. 1.1024 4.7094][0.0000 0.9055 4.9063][0.0000 -0.0945 

4.9063][0.0000 -0.1112 4.8993][0.0000 -0.1181 4.8826][0.0000 
0.0000 3.9138][72.0000 0.0000 3.9138][0.0000 0.0000 
0.0000][72.0000 0.0000 0.0000][0.0000 -20.0000 0.0000][72.0000 
_-20.0000 0.0000][0.0000 20.0000 0.0000][72.0000 20.0000 0.0000] 


To build the cireular flange stiffened plates, ensure a mesh is ereated for the curve 
at the top of the web with topology BAR2. The curve will serve as the flange once the 
appropriate beam section is created. Using the properties menu, create and name a ID 
beam property set. Use the general section (CBEAM) option. Use the input properties 
button and then the create sections beam library button to create and emplace the 
appropriate flange onto the web. Use the beam library to create and name a flange by 
selecting the standard shape cross-section using the solid circular shape button with the 
appropriate radius dimension. In the general beam (CBEAM) input properties 
subordinate menu, enter the appropriate material name, section name (just created in the 
create sections beam library), and bar orientation. 

Table 48 Data for circular flange stiffened plate model Ic (units in inches) 

See Table 7 model [0.0000 0.0000 3.730197][72.0000 0.0000 3.730197][0.0000 
Ic for parameters. 0.0000 0.0000][72.0000 0.0000 0.0000][0.0000 -20.0000 
0.0000][72.0000 -20.0000 0.0000][0.0000 20.0000 
_ 0.00001[72.0000 20.0000 0.0000] _ 
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Table 49 Data for circular flange stiffened plate model 2c (units in inches) 


See Table 8 model 
2c for parameters. 

[0.0000 0.0000 3.768468][72.0000 0.0000 3.768468][0.0000 

0.0000 0.0000][72.0000 0.0000 0.0000][0.0000 -20.0000 
0.0000][72.0000 -20.0000 0.0000][0.0000 20.0000 

0.0000][72.0000 20.0000 0.0000] 

Table 50 Data for circular flange stiffened plate model 3c (units in inches) 

See Table 9 model 
3c for parameters. 

[0.0000 0.0000 3.831455][72.0000 0.0000 3.831455][0.0000 

0.0000 0.0000][72.0000 0.0000 0.0000][0.0000 -20.0000 
0.0000][72.0000 -20.0000 0.0000][0.0000 20.0000 

0.0000][72.0000 20.0000 0.0000] 

Table 51 Data for circular flange stiffened plate model 4c (units in inches) 

See Table 10 
model 4c for 
parameters. 

[0.0000 0.0000 4.320748][72.0000 0.0000 4.320748][0.0000 

0.0000 0.0000][72.0000 0.0000 0.0000][0.0000 -20.0000 
0.0000][72.0000 -20.0000 0.0000][0.0000 20.0000 

0.0000][72.0000 20.0000 0.0000] 

Table 52 Data for circular flange stiffened plate model 5c (units in inches) 

See Table 11 
model 5 c for 
parameters. 

[0.0000 0.0000 4.359019][72.0000 0.0000 4.359019][0.0000 

0.0000 0.0000][72.0000 0.0000 0.0000][0.0000 -20.0000 
0.0000][72.0000 -20.0000 0.0000][0.0000 20.0000 

0.0000][72.0000 20.0000 0.0000] 

Table 53 Data for circular flange stiffened plate model 6c (units in inches) 

See Table 12 
model 6c for 
parameters. 

[0.0000 0.0000 4.422005][72.0000 0.0000 4.422005][0.0000 

0.0000 0.0000][72.0000 0.0000 0.0000][0.0000 -20.0000 
0.0000][72.0000 -20.0000 0.0000][0.0000 20.0000 

0.0000][72.0000 20.0000 0.0000] 


To build the T-flange stiffened plate models, the same data in Table 54 is used to 
build T-flange stiffened plate models It, 2t, and 3t. The same data in Table 55 is used to 
build T-flange stiffened plate models 4t, 5t, and 6t. Ensure a mesh is ereated for the 
curve at the top of the web with topology BAR2. The curve will serve as the flange once 
the appropriate beam section is created. Using the properties menu, create and name a 
ID beam property set. Use the general section (CBEAM) option. Use the input 
properties button and then the create sections beam library button to create and emplace 
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the appropriate flange onto the web. In the beam library, ereate and name a standard 
shape eross-section using the solid rectangular shape button with the appropriate 
dimensions. In the general beam (CBEAM) input properties subordinate menu, enter the 
appropriate material name, section name (just created in the create sections beam library), 
and bar orientation (start by using the vector (0,1,0) and adjust with (l,0,0) or (0,0,l) 
as needed). 


Table 54 Data for T-flange stiffened plate models It, 2t, and 3t (units in inches) 


See Table 7 model 
It, Table 8 model 

2t, or Table 9 
model 3t as 
appropriate for 
parameters. 

[0.0000 0.0000 3.4414][72.0000 0.0000 3.4414][0.0000 0.0000 
0.0000][72.0000 0.0000 0.0000][0.0000 -20.0000 0.0000][72.0000 
-20.0000 0.0000][0.0000 20.0000 0.0000][72.0000 20.0000 0.0000] 

Table 55 Data for T-flange stiffened plate model 4t, 5t, and 6t (units in inches) 

See Table 10 
model 4t, Table 11 
model 5t, or Table 

12 model 6t as 
appropriate for 
parameters. 

[0.0000 0.0000 4.031951][72.0000 0.0000 4.031951][0.0000 

0.0000 0.0000][72.0000 0.0000 0.0000][0.0000 -20.0000 
0.0000][72.0000 -20.0000 0.0000][0.0000 20.0000 

0.0000][72.0000 20.0000 0.0000] 


97 





THIS PAGE INTENTIONALLY LEET BLANK 


98 



LIST OF REFERENCES 


' Bleich, F., Buckling Strength of Metal Structures, McGraw-Hill Book Co., 1952. 

2- Chajes, Alexander, Principles of Structural Stability Theory, Prentice-Hall, Inc., 
1974. 

Chou, S. K. G., Design of Longitudinally Stiffened Plating in Compression with 
Particular Reference to Torsional Buckling, Ph.D. Thesis, Civil Engineering Dept., 
Imperial College, London, August 1997. 

4- Corus Group, “Bulb flats”, [http://www.oorus-specialprofdes.com/pdf/Bulbflat.pdf]. 
September 2002. 

5- Danielson, D. A., “Theory of Shell Stability”, Thin-Shell Structures, pp. 45-58, 1974. 

6 Danielson, D. A., Kihl, D.P., and Hodges, D.H., “Tripping of Thin-Walled Plating 

Stiffeners in Axial Compression”, Thin-Walled Structures, Vol. 10, pp. 121-142, 
1990. 

Danielson, D. A., Crioelli, A.S., Frenzen, C.L., and Vasudevan, N., “Buckling of 
Stiffened Plates Under Axial Compression and Lateral Pressure”, Int. J. Solids 
Structures, Vol. 30, No. 4, pp 545-551, 1993. 

^ Danielson, D. A., Steele, C. R., Fakhroo, F. and Cricelli, A. S., Stresses in Ship 
Plating, Naval Postgraduate School Technical Report NPS-94-008, 1994. 

^ Danielson, D.A., Analytieal Tripping Loads for Stiffened Plates, Int. J. Solids 
Structures, Vol. 32, No. 8/9, pp. 1319-1328, 1995. 

Donaldson, B. K., Analysis of Aircraft Structures: An Introduction, McGraw-Hill, 
Inc., 1993. 

Henrici, P., Applied and Computational Complex Analysis, Volume I, Power Series- 
Integration-Conformal Mapping-Location of Zeros, John Wiley & Sons, 1974. 

'2- Jaeger, L. G., Elementary Theory of Elastic Plates, Pergamon Press, 1964. 

'3- Langhaar, H. L., Energy Methods in Applied Mechanics, John Wiley and Sons, Inc., 
1962. 

'4 Naval Postgraduate Sehool Teehnical Report NPS-MA-96-002, Buckling of Ship 

Grillages, by D. Danielson and D.P. Kihl, September 1996. 

99 



Naval Postgraduate School Technical Report NPS-MA-97-005, Buckling of Ship 
Grillages - Part II, by D. Danielson and D.P. Kihl, September 1997. 

Oden, J. T., Mechanics of Elastic Structures, McGraw-Hill Book Co., 1967. 

Paik, J. K. and Kim, B. J., “Ultimate Strength Formulations for Stiffened Plates 
Under Combined Axial Load, In-Plane Bending and Lateral Pressure: A Benchmark 
Study”, Thin-Walled Structures, Vol. 40, pp.45-83, 2002. 

Pilkey, W. D. and Chang, P.Y., Modern Formulas for Static and Dynamics: A Stress- 
and-Strain Approach, McGraw-Hill Book Co. 1978. 

Sokolnikoff, 1. S., Mathematical Theory of Elasticity, 2"‘* ed., McGraw-Hill Book Co., 
Inc., 1956. 

20- Stakgold, L, Green’s Functions and Boundary Value Problems, 2“‘* ed., John Wiley & 
Sons, Inc., 1998. 

21- Timoshenko, S. P. and Gere, J. M., Theory of Elastic Stability, 2“‘* ed., McGraw-Hill 
Book Co., 1961. 

22- Timoshenko, S. P. and Goodier, J. N., Theory of Elasticity, 3’^‘* ed., McGraw-Hill 
Book Co., 1970. 

22- Ugural, A. C., Stresses in Plates and Shells, McGraw-Hill, Inc., 1981. 

24- Washizu, K., Variational Methods in Elasticity and Plasticity, 3'^‘* ed., Pergamon 
Press, 1982. 


100 



INITIAL DISTRIBUTION LIST 


1. Defense Teehnieal Information Center 
Ft. Belvoir, Virginia 

2. Dudley Knox Library 
Naval Postgraduate Sehool 
Monterey, California 

3. Dr. Don Danielson 

Naval Postgraduate Sehool 
Monterey, California 

4. Dr. Young Kwon 

Naval Postgraduate Sehool 
Monterey, California 

5. Dr. David Canright 

Naval Postgraduate Sehool 
Monterey, California 

6. Dr. Fariba Fahroo 

Naval Postgraduate Sehool 
Monterey, California 

7. Dr. Wei Kang 

Naval Postgraduate Sehool 
Monterey, California 

8. Dr. David Kihl 

Naval Surfaee Warfare Center 
West Bethesda, Maryland 

9. Geoff Lofthouse 
Corns Group 

Skinningrove, United Kingdom 


101 



